Unclassified 

SECURITY  CLASSIFICATION  of  THIS  PACE  (Whmn  Dmtm  Enfrrd) 


J ONRijCR212-224-3 


lc 


4.  TITLE  (end  Subtitle) 

SELF  OPTIMIZING  FLEXIBLE  TECHNOLOGY  \yiNG 
PROGRAM,  Semispan  Wind  Tunnel  Test  of  a General- 
]_.ized  Geometry  Wing  Panel 

Model  Design  and  Aerodynamic 
Analysis  , 


REPORT  DOCUMENTATION  PAGE 


3RT  NUMBER 


W 


2.  GOVT  ACCESSION  NO 


7.  AUTMORf.) 

S./Levinsky 
R.  H./Schappelle 
A.  A .1  McClain 

a.  PERFORMING  OHU'JInU  ATION  N AME  AND  ADDRESS 

General  Dynamics  Corporation  Convair  Division 
5001  Kearny  Villa  Road, 

San  Diego,  California  92112 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


Office  of  Naval  Research  - Department  of  the  Navy 
800  North  Quincy  St. 

Arlington,  VA  22217 


I*  MONITORING  AGENCY  NAME  6 ADORESSOI  ditto, rnn!  from  Controlling  Oltico) 


j ^ F/r  / ' / U 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3-  RECIPIENT'S  CATALOG  NUMBE. 


5 , TYPE  OP  REPORT  ft  PERIOO  COMB  HEP 

^ Interim  Technical  Repwrt. 
15  Apr#  M76  — 14  Feb,  *177, 


ASD-NSC-77-001^j 

15  j— 


“v* 

avnuMT-rwcBEHi.r 


USER 


^N00014-76-C  -0742  J 


10  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A WORK  UNIT  NUMBERS 

02241N 

RF  41-411-801 

NR  212-224 


12  REPORT  DATE 

February  14,  1977 


13.  NUMBER  OF  PAGES 


15.  SECURITY  CLASS  tot  thia  report) 

Unclassified 


IS*.  OECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 


16  DISTRIBUTION  STATEMENT  (ol  thlo  Rmporl) 


(jjj)  /¥  V&W  7/ 

Approved  for  Public  Release;  Distribution  Unlimited  ^ ^ ^ | 


17.  DISTRIBUTION  STATEMENT  (ol  the  abatrmet  entered  In  Block  20,  II  dlllerent  Irom  Report) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  reveree  aide  II  neceaaery  end  Identity  by  block  number) 

Optimization  Mission  Adaptive  Wing 

Wind  Tunnel  Models  Nonlinear  Programming 

Wind  Tunnel  Testing  Transonic  Flow 

Variable  Camber  Wing  Aerodynamics 


Interactive  Testing 
Computer  Controlled  Testing 


ABSTRACT  (Contlnn,  nn  rmvtrra,  aid,  II  naraaaary  and  Identity  hy  block  numbar} 

Various  aspects  of  the  design  and  stress,  flutter,  aerodynamic  and  computer  analyses 
are  presented  for  a l/6^h  scale  semispan  Self-Optimizing  Flexible  Technology  (SOFT) 
Wing  Wind  Tunnel  Model  representative  of  the  AFTI-111  configuration,  to  be  tested 
at  AEDC  tunnel  16T.  The  model  features  twelve  independent  hydraulic  actuator 
systems  which  are  designed  to  conformally  vary  leading  edge  radius,  leading  edge 
deflection  about  the  15%  and  25 % chord  lines,  upper  surface  humping,  and  trailing 


dd 


FORM 
AN  7] 


1473 


EDITION  OF  I NOV  6S  IS  OBSOLETE 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whmn  Dmtm  Enl»r»«J) 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE(W7»«n  Dmtm  EnCfd) 

'"edge  deflection  about  the  65%  and  80%  chord  lines  at  an  inboard  and  outboard  station. 
The  model  is  designed  to  be  tested  interactively  under  computer  control  in  such  a 
manner  so  as  to  optimize  any  one  of  eighteen  possible  aerodynamic  (or  structux-al) 
functions,  while  satisfying  equality  (or  inequality)  constraints  with  the  remaining 
seventeen.  These  functions  are  to  be  measured  "on  line"  by  either  the  wind  tunnel 
balance,  potentiometers,  wing  pressure  orifices  or  strain  gages 

Simulations  of  the  wing  optimization  procedure  are  presented  using  a highly 
simplified  set  of  nonlinear  algebraic  equations  to  generate  the  aerodynamic  functions. 
The  simulations  show  that  convergence  may  be  obtained  for  various  optimization 
problems,  even  when  randomly  generated  testing  type  errors  are  included,  provided 
however,  that  the  various  parameters  in  the  optimization  computer  program  have  been 
properly  tuned  for  the  particular  optimization  problem.^ 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfTFhen  f)»r«  F.ntr'ed) 


hr 


J 


FOREWORD 


This  research  program  was  undertaken  by  the  Convair  Division  of  General 
Dynamics  Corporation,  San  Diego,  California,  under  Office  of  Naval 
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NOMENCLATURE 


Ajj  actuator  position  (counts) 

A1  aspect  ratio 

a 2D  lift  curve  slope 

b full  span  of  model 

bj  effective  semispan  of  i th  wing  station 

Cq  drag  coefficient 

CDp  friction  drag  coefficient  of  entire  configuration 

C^ind  Educed  drag  coefficient 

AC f>p  additional  profile  drag  coefficient 

ACdri  additional  compressibility  drag  coefficient 
Cl  lift  coefficient 

Cl  lift  coefficient  corrected  for  tail  lift  required  to  trim 
C;vi  pitching  moment  coefficient 

Cp  pressure  coefficient 

Crb  ro°t  bending  moment  coefficient  (=  Mbb/i813) 

CSC  creepless  supercritical  airfoil  section 

c local  wing  chord  (measured  streamwise  with  A - 16°, 

also  sound  speed 

Cj  average  chord  of  i' th  actuator 

c mean  aerodynamic  chord 

e aerodynamic  efficiency  factor 

F Lagrangian  function  being  minimized 

g vector  of  active  constraints 

H Hessian  matrix 

H (x)  unit  step  function 

step  size  factor  used  during  simultaneous  mode 
A Kj  partial  span  lift  factors 

tail  arm 

M free  stream  Mach  number  (unless  subscripted  otherwise), 

also  moment 


ix 


NOMENCLATURE  (Cont’d) 


MS  model  scale 

Pj  . compressibility  drag  correction  factor 

q free  stream  dynamic  pressure  (unless  subscripted  otherwise) 

R restoration  criterion  function 

Rs  i leading  edge  radius  of  curvature 

Re  Reynolds  number  based  on  c 

S wing  reference  half  area  (unless  subscripted  otherwise); 

also  vector  of  search  directions 

Sij,  Sq.  search  direction  components  for  minimizing  F 

Sr  search  direction  vector  for  minimizing  R 

t airfoil  thickness  or  maximum  thickness 

VV52  TACT  wing  52 

x,  z axial  and  vertical  positions 

x , axial  position  of  centroid  of  lift  of  fuselage  and  glove 
GLOVE 

o angle  of  attack  (degs) 

tail  incidence  angle  (degs) 
y ratio  of  specific  heats  (y=  1.4) 

A incremental  value 

6. j actuator  deflection  angles  (degs) 

et  tail  downwash  angle  (degs) 

T|jj  actuator  aerodynamic  effectiveness  factors 

0- j , Qa  nondimensional  transformed  coordinates  used  in 
optimization  program 

fij  nondimensional  actuator  chords 

j\  sweep  angle  (degs) 

A vector  of  Lagrange  multipliers  used  in 

optimization  program 

a stress  level 
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NOMENCLATURE  (Cont'd) 

0 objective  (or  merit)  function  being  minimized 

potentiometer  resistance  (ohms) 

V gradient  operator 

SUBSCRIPTS 

c value  beyond  stall  at  which  lift  coefficient  stops 

decreasing 

i,  j indicies  referring  to  actuator  system 
i = 1 (inboard),  2 (outboard) 
j = 1 (L.E.  radius),  2 (deflection  about  15%c), 

3 (deflection  about  25%c),  4 (upper  surface  humping), 
5 (deflection  about  G5%c),  6 (deflection  about  80%  c) 
k test  point  index  during  simultaneous  mode 

L lower  wing  surface 

MAX  maximum 

MC  moment  center 

MIN  minimum 

P location  of  wing  pivot,  also  pressure 

RB  root  bending 

r restoration  direction 

S.F.  safety  factor 

STALL  value  at  stall 

TEj  trailing  edge  value  at  station  i 
TWj  wing  twist  at  station  i 

t tail 

U upper  wing  surface 

OLj  zero  lift  value  for  station  i 

a 9/S  a 
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SUSCRIPTS 


6 3/9  6 

00  free  stream 

SUPERSCRIPTS 
T transposed  matrix 

i * nominal  value,  also  best  value 

See  Table  5. 1 for  common  variables  in  the  optimization 
computer  program. 

1 


xii 


The  Self  Optimizing  Flexible  Technology  (SOFT)  Wing  program  deals  with  the  design, 
fabrication  and  testing  of  a SOFT  Wing  wind  tunnel  model  and  associated  computer 
and  control  systems  at  AEDC  Tunnel  1GT.  The  program  is  being  conducted  by 
General  Dynamics  Convair  Division  (GDC)  and  the  Arnold  Engineering  and  Development 
Center  (AEDC)/Aerodynamics  Research  Organization  Inc.  (ARO)  on  a cooperative 
basis.  Funding  is  being  provided  by  U.S.  Navy  (Office  of  Naval  Research),  and  the 
U.S.  .Air  Force  (Flight  Dynamics  Laboratory  and  AEDC/ARO). 

1. 1 PROGRAM  STATUS 

In  order  to  minimize  the  overall  technical  and  financial  risks  and  to  remain  within 
budgetary  limits,  the  program  has  been  incrementally  funded  and  has  been  divided 
into  two  distinct  phases: 

(I)  Model  design,  stress,  flutter  and  aerodynamic  analyses,  computer 
programming  and  test  simulation  ~ GDC 

(II)  Model  fabrication,  calibration,  test  planning,  test  support  and 
balance  calibration  ~ GDC 

and 

Control  system,  instrumentation,  installation, 
computer  program  interface  and  testing  - ARO 


Specific  tasks  were  to  have  been  carried  out  during  each  phase  of  the  program 
according  to  the  schedule  given  in  the  original  technical  proposal  (Ref.  1).  The 
schedule  has  been  subsequently  revised  as  shown  in  Table  1.1.  The  decision  to  go- 
ahead  with  Phase  (II)  was  made  at  the  prefabrication  conference,  held  at  WPAFB, 
Ohio,  on  14  October  1976,  (Ref.  2).  This  meeting  is  labeled  as  meeting  #4  under 
Task  4.  2.  8b  in  Table  1. 1. 

The  status  of  the  various  Phase  (II)  tasks  is  currently  as  shown  in  Table  1.2. 
All  Phase  (II)  tasks  are  either  fully  funded,  or  in  the  final  stages  of  negotiation. 


except  tor  Tasks  4.  2.  6d  and  4.  2.  8d  dealing'  with  the  aerodynamic  evaluation  of  the 
test  data  and  the  preparation  and  presentation  of  a final  report.  It  is  planned  to  submit 
a proposal  for  these  tasks  prior  to  the  wind  tunnel  test,  now  set  for  July  14,  1077  as 
shown  in  Table  1.  1.  The  pretest  conference  'meeting  #5  under  Task  4.  2.  8b)  was 
held  at  GDC.  San  Diego,  Calif,  over  the  period  13-16  Dec.  l‘>76  Ref.  3).  Details 
ot  the  various  Phase  (II)  activities  may  be  obtained  from  the  pretest  report  already 
submitted  under  Task  4.  2.6  c (Ref.  4). 

The  present  interim  technical  report  deals  only  with  the  Phase  I effort, 
which  was  sponsored  by  the  U.S.  Navy  (ONR).  The  report  is  scheduled  under 
Task  4. 1.7  (iii)  of  Table  1.  1.  Loads,  stress  and  flutter  data  and  final  design 
drawings  are  being  delivered  under  separate  cover  (Refs.  5-7),  in  accordance  with 
Task  4. 1.7  (ii)  of  Table  1. 1.  Tapes  and  listings  of  the  optimization  computer 
program  developed  under  Task  4. 1.5  (i)  were  submitted  previously  to  ARO  during  the 
pretest  conference  (Ref.  3). 


1.  2 PROGRAM  DESCRIPTION  AND  BACKGROUND 

The  General  Dynamics  Convair  SOFT  Wing  concept  is  a wing  winch  adapts  its  shape 
conformally  (continuous  chordwise  and  spanwise  variations  of  camber,  thickness  and 
twist)  in  a manner  to  maximize  or  minimize  various  merit  functions  (e.g. , maximum 
load  factor,  minimum  drag,  maximum  lift/drag  ratio),  while  satisfying  constraints 
(e.g.  , constant  lift,  trim,  the  avoidance  of  excessive  buffet  or  flow  separation,  the 
avoidance  of  excessive  structural  stresses).  The  SOFT  Wing  thereby  differs  from 
other  variable  camber  systems  in  that  it  not  only  employs  mechanisms  to  conformally 
vary  the  geometry,  but  also  automatically  adapts  its  shape  to  each  flight  (or  test) 
condition.  Ideally,  the  wing  shapes  are  optimum  (within  the  limitations  of  the 
structural  flexibility  and  articulation)  and  include  all  of  the  various  viscous  and 
in  viscid  flow  effects. 

Although  the  eventual  application  of  the  concept  may  be  to  real  time  flight 
conditions,  the  present  program  deals  with  the  wind  tunnel  testing  and  optimization 
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of  a l/6'th  scale  semi-span  SOFT  Wing  model  representative  of  the  AFTI-111 
configuration.  The  test  has  been  scheduled  in  AEDC  tunnel  1GT  during  the  summer 
of  1977,  and  a test  period  of  128  hours  occupancy  is  planned.  Purposes  of  the 
investigation  are  as  follows: 

(i)  to  verify  the  structural  integrity,  flexibility  and  controllability 
of  the  model  under  transonic  flow  conditions, 

(ii)  to  determine  the  suitability  of  the  closed-loop  computer  control  system, 

(iii)  to  verify  the  convergence  and  stability  characteristics  of  the  optimiza- 
tion computer  algorithms  under  conditions  of  transonic  flow  and  with 
up  to  thirteen  different  control  variables, 

(iv)  to  arrive  at  optimum  wing  shapes  for  a given  test  condition,  merit 
function  and  set  of  constraints,  utilizing  the  maximum  articulation 
capabilities  of  the  model. 

(v)  to  evaluate  possible  improvements  in  performance,  maneuverability 
and  load  control  which  might  accrue  to  the  AFTI-111  aircraft  through 
the  utilization  of  a mission  adaptive  wing  with  either  full  model 
articulation  or  more  limited  articulation  dictated  by  current  AFTI-111 
designs. 

(vi)  to  demonstrate  and  use  the  General  Dynamics  Fort  Worth  (GDFW) 
stereophotographic  system  for  measuring  wing  shape  under  airload 
conditions  in  the  tunnel. 

Recent  test  results  obtained  by  GDFW  in  connection  with  the  AFTI-111 
variable  camber,  mission  adaptive  wing  program,  involving  a wind  tunnel  model 
with  interchangeable  parts,  have  shown  that  the  optimum  leading  edge  and  trailing 
edge  camber  and  deflection  are  interdependent  and  sensitive  to  flight  condition 
and  static  margin,  (Ref.  8).  The  AFTI-111  tests  seem  to  show  that  a large  saving 
in  wind  tunnel  test  time  might  well  have  been  achieved  through  use  of  a SOFT 
Wing  type  model.  This  would  be  especially  true  as  the  number  of  degrees  of 
freedom  and  extent  of  articulation  increase,  and  because  the  precision  of  control 
for  the  SOFT  Wing  permits  essentially  a continuous  spectrum  of  settings. 
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Under  a previous  contract  for  the  U.S.  Navy  (ONR),  GDC  has  succeeded  in 
fabricating  and  testing  a two-dimensional  flexible  airfoil  model  which  could  be 
optimized  in  the  tunnel  in  a manner  similar  to  that  planned  for  the  present  SOFT 
Wing  model.  The  2D  model  was  tested  under  both  low-speed^)  and  transonic  ^ 
conditions , and  employed  five  pair  of  hydraulic  actuators  to  independently  vary 
leading  edge  radius,  leading  edge  incidence,  lower  surface  humping,  lower  surface 
aft  camber  and  trailing  edge  deflection.  During  the  transonic  tests  (M=0.  75)  a link- 
up of  analog  and  digital  computers  was  used.  The  airfoil  shape  and  angle  of  attack  were 
controlled  "on  line"  to  minimize  drag,  subject  to  constraints  on  lift  and  trailing  edge 
pressure,  using  the  gradient  projection  optimization  algorithm.  Satisfactory 
convergence  was  found,  even  with  flow  separation;  however,  because  of  the  inter- 
mittent blow-down  type  operation  of  the  transonic  test  facility,  a sequence  of 
individual  inns  was  required  to  complete  a single  optimization  problem. 

The  new  3D  SOFT  Wing  model  employs  twelve  independent  hydraulic  actuator 
systems  of  similar  design  to  those  for  the  two-dimension  model.  The  actuators  are 
designed  to  perform  the  following  variations: 

A (nose  radius,  inboard  station) 

A (nose  deflection  about  15%  chord  line,  inboard  station) 

A (nose  deflection  about  25%  chord  line,  inboard  station) 

it) 

A j ( (upper  surface  humping,  inboard  station) 

A}-  (trailing  edge  deflection  about  65%  chord  line,  inboard  station) 

A (trailing  edge  deflection  about  80%  chord  line,  inboard  station) 

16 

A - A,„  (same  as  A,,  - A . except  at  outboard  stations). 

21  26  11  16 

It  is  planned  to  employ  similar  computer  controlled  testing  and  optimization  techniques 
with  the  new  model  as  used  for  the  2-D  investigations.  Hence  the  new  actuator 
systems,  position  readout  instrumentation,  control  valves  and  computer  link  ups 
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are  similar  to  those  for  the  2D  model,  except  that  the  number  of  control  channels  has 
been  increased  to  handle  the  larger  number  of  independent  control  variables  (13 
including  angle  of  attack).  Other  differences  relate  to  the  increased  volume  and 
pressure  of  the  hydraulic  system,  to  the  incorporation  of  hydraulic  shut-off  valves  to 
lock  the  actuator  positions  in  the  event  of  a failure  of  some  component,  and  to  the 
increased  accuracy  in  positioning  the  actuators  afforded  by  the  new  DAC  and  ADC 
components.  The  design,  assembly  and  demonstration  of  the  new  control  system 
(Task  4.  2.  2 of  Table  1. 1)  is  being  performed  by  AEDC/ARO  under  Phase  II,  and 
will  not  be  reviewed  further  herein,  except  to  note  that  the  2D  model  is  being  used 
by  ARO  to  check  out  the  new  control  system  prior  to  completion  of  the  3D  model. 

The  forthcoming  tests  will  use  the  computer  elements  of  the  AEDC  closed- 
loop  control  system,  which  was  developed  in  the  Propulsive  Wind  Tunnel  (PWT) 
Facility,  and  is  described  in  References  12  and  13.  Because  of  the  continuous 
operation  of  tunnel  16T,  it  is  planned  to  exercise  continuous  computer  control  of 
the  model  during  wing  optimization.  Hence,  under  ideal  circumstances,  a single 
run  in  the  wind  tunnel  would  be  sufficient  to  arrive  at  an  optimum  wing  shape  for  a 
prescribed  flight  condition  (Mach  number  and  Reynolds  number)  and  given  merit 
function  and  constraints. 

Because  of  the  increased  number  of  independent  variables,  and  because 
of  the  greatly  increased  number  and  varie'  of  the  constraints  for  the  3D  wing 
optimization  problem,  with  respect  to  that  for  the  2D  model,  the  optimization 
prog-am  allows  the  merit  function  to  be  optimized  to  be  any  one  of  eighteen  functions 
(or  combinations  of  functions)  which  are  measured  by  either  the  wind  tunnel  balance, 
the  model  potentiometers,  pressure  orifices,  or  strain  gages.  Any,  or  all  of  the 
remaining  seventeen  functions  may  be  specified  as  equality  or  inequality  constraints. 
Among  the  new  constraints  which  have  been  used  in  the  program  are  inequality 
constraints  which  limit  the  maximum  differential  deflection  between  corresponding 
inboard  and  outboard  systems  (termed  differential  twist  constraints)  and  inequality 
constraints  on  maximum  rolling  moment  coefficient  (for  load  control) 
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and  on  maximum  wing  root  bending  stress.  In  addition,  equality  constraints  on 
pitching  moment  coefficient  may  be  introduced  via  table  lookup.  This  permits  the 
optimization  problem  to  account  for  effects  of  trim,  without  exceeding  the  maximum 
horizontal  tail  effectiveness,  even  through  the  tail  incidence  on  the  model  is  held 
fixed  during  the  optimization.  The  new  optimization  problem  allows  up  to  fourteen 
independent  variables  consisting  of  the  twelve  actuators  A-,  i = 1,  2;  j = 1,  2,  . . .fi, 
the  tail  incidence  angle  and  angle  of  attack  a.  Any  one  or  several  of  the 
independent  variables  may  be  held  fixed  during  an  optimization  problem,  in  order  to 
simulate  variable  camber  systems  with  reduced  articulation,  such  as  presently 
envisioned  for  the  AFT1-111  (Ref.  8). 

The  new  optimization  program  has  been  tested  extensively  for  convergence 
in  a series  of  simulated  optimization  problems.  The  aerodynamic  characteristics 
used  in  the  program  were  obtained  from  a set  of  nonlinear  algebraic  equations,  termed 
"Pseudo  Aerodynamic  Module"  (PAM),  which  were  linked  to  the  optimization  program 
for  the  simulations.  PAM  provides  aerodynamic  coefficients  as  functions  of  the 
14  independent  variables,  Mach  number  and  Reynolds  number.  Emphasis  in  the  PAM 
was  placed  on  speed  of  computation  and  on  a realistic  representation  of  the  non  linearities 
expected  in  the  test  data,  rather  than  on  accuracy,  in  order  to  hold  computational  time 
in  the  simulations  to  an  absolute  minimum  Testing  type  errors  and  inaccuracies  were 
added  to  PAM  to  further  improve  the  realism  of  the  simulations. 

A summary  of  the  design,  including  the  support  system,  fuselage,  leading 
edge  and  trailing  edge  actuation  systems  and  instrumentation  lias  been  included  in 
the  next  section.  A brief  summary  of  the  stress,  loads,  and  flutter  considerations 
is  given  in  Section  :i,  whereas  the  evaluation  of  a plastic,  kinematic  model  of  the 
SOFT  Wing  is  found  in  Section  4.  The  mathematical  basis  for  the  optimization  program, 
the  defining  equations  in  PAM,  and  results  of  test  simulation  computations  carried 
out  on  the  DEC  10  computer  are  presented  in  Section  5 Various  aerodynamic 
computations  already  carried  out  and  planned  in  the  future  for  comparing  the  aero- 
dynamic characteristics  of  various  desired  (target,  AFTI-111  and  FACT  1 airfoils 


with  the  closest  approximating  (actual)  SOFT  Wing  airfoils  are  found  in  Section  G. 
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AHI.f.  1 1 OVERALL  SCHEDU LE  SOFT  WING  PROGRAM,  PHASES  I AND  n REVISED 
1)l:('  31.  FOR  I ACT  TEST  AT  •VUES,  JULY  i5  A ~UC  TEST  DATE 


Phase 

and 

Contract 

Designation* 


Basic 

ONF 

Desigr.  i 
Stud> 


Months  From  Go-  Ahead  1 

6 7b  0 10  11  12  13  14  15  13  17 


F M A I M J 


II  4.2.2 
4.2.3 


AEDC  4.2.6 


Model  Predesign 
Model  Design 
Stress/Loads/  Flutter 
Kinematic  Model 

(i)  Fabrication 
(U)  Demonstration 
Computer  Simulation  Studies 

(1)  Expand  Optimization  Program 

(ii)  Simplified  Aero  Module 

(iii)  Conduct  Test  Simulation 
Aerodynamic  Analyses 

(i)  Compare  with  F-lll  TACT 

(ii)  Evaluate  Actuator  Effectiveness 
Documentation 

(i)  Progress  Letters 

(ID  Loads/Stress  Flutter  Report 

(ill)  Interim  Technical  Report 

Model  Design  Modification  for  TACT 

Program 

(!)  • (tv)  Design 

iv)  Computer  Programming 

liaison  and  Prefabrication  Activities 

Purchase  of  Long  Lead  Time  Materials 

Model  Fabrication  and  Assembly 

а.  Wing  Leading  Edge 

б.  Wing  Trailing  Edge 

c.  Materials 

d.  Calibration  Tooling 
Shipment  to  AEDC 

e.  Fuselage  Location  

Model  Control  System 

Model  Calibration 

a.  Model  Setup  at  AEDC 

b.  Hydraulics  t Instrumentation 

c.  Skin -off  Calibration 
d Sldn-on  Calibration 
e.  Stereo -Photography 
i.  Install  Fuselage 
Computer  Programming 

a.  Data  Reduction  Program 

b.  Optimization  Program 

c.  Simplified  Aero  Module 

d.  Test  Simulation 

e.  Final  Link-up  and  Demonstration 
Test  Planning  and  Operation 

a.  Pretest  Planning  & Coordination 

b.  Test  Support 

Data  Evaluation  and  Documentation 

a.  On-Line  Data  Analysis* 

b.  Monthly  Progress  Letters 

c.  Pretest  Report 

d.  Data  Eva  1/ Draft  Technical  Report* 

e.  Software  Listing 

f.  Revised  Technical  Report* 

Wind  Tunnel  Balance  Calibration 

a.  Pretest  Callb.  at  GDC.  S.D. . Calif. 

b.  Check  Calibrate  at  AEDC 

c.  Post- test  Callb.  at  GDC.  S.D. . Calif. 
Milestones  L Meetings 


U 


n 

I L II 
II.  5 
Unfunded 

a.  Start  Date 

b.  Meetings** 

c.  Trips  lo  AEDC  <Flw  Total1 

d.  Presentation  of  Final  Report* 
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4.2.9  AEDC  Supplemtait-il  Fabrication  Tasks 

D MOD 

a.  Final  Assembly  L Functional  Checkout 

AEDC/ 

b Fuselage  & Glove  Rework 

FDL 

c Materials  fc  Instrumentation 

? 

d.  Engineering  Liaison 

12b  Hour  Test 
AEDC  16T 


SolM  Bar.  GDC  Task  for  U.S.  Navy  (ONR) 

Open  Bar  GDC  Assist  Task  for  AEDC  (ARO  Inc.) 

Cross  Hatched  Bar  Task  to  be  jxrformed  solely  by  AEDC  (ARO  Inc.) 
• Aspects  dealing  with  post-test  data  evaluation  and  final  report 
are  presently  unfunded 

••  Meeting?*!  - *5  are  described  tn  Section  2.0  under  task  4.2. 8.  b 


2.0  SUMMARY  OF  MODEL  DESIGN 


The  wind  tunnel  model  is  a l/fi'th  scale  semi-span  version  of  the  F-lll  TACT 
configuration  with  a SOFT  Wing  capability  employing  twelve  remotely  controllable 
hydraulic  actuator  systems,  separable  front  and  rear  spars,  and  three  flexible 
steel  skin  panels  with  sliding  joints.  The  spar  box  and  skin  sections  have  been 
formed  to  correspond  to  the  TACT  W52  supercritical  airfoil  shapes  (unloaded)  at 
two  control  stations  (0.348  b/2  and  0.913  b/2)  for  a given  set  of  actuator  positions, 
to  be  determined  as  part  of  the  calibration.  The  wing  and  fuselage  are  designed  for 
testing  at  sweep  angles  of  26°  and  58°,  with  or  without  a horizontal  tail  at  several 
fixed  incidence  angles.  Additional  sweep  angles,  viz. , ICC,  450  and  G5°  may  be  set 
by  making  minor  modifications  to  the  fuselage  glove  section.  However,  trailing 
edge  clearance  problems  will  limit  the  maximum  trailing  edge  deflection  of  the 
inboard  actuators  at  sweep  angles  at  or  beyond  45°. 

2.1  SUPPORT  SYSTEM 

The  3-D  SOFT  Wing  model  will  be  mounted  horizontally  from  the  east  sidewall 
and  will  be  installed  as  depicted  in  Figure  2. 1.  The  model  will  be  supported  by 
the  NASA  Ames  C-120-16.36  wind  tunnel  balance,  to  be  procured  by  the  Air  Force 
FDD  and  modified  and  calibrated  by  GDC  during  the  Phase  II  effort.  The  fuselage 
(shown  in  Figure  2. 1)  is  metric  and  is  mounted  off  a nonmetric  image  plate,  which 
rotates  with  the  fuselage.  The  gap,  of  approximately  1/4  inch,  between  the  fuselage 
and  the  image  plate  shall  be  sealed  with  a fixed  flexible  seal  designed  to  minimize 
any  transmission  of  loads.  The  image  plate  is  mounted  off  of  a reflection  plane  which 
extends  out  approximately  thirty  inches  from  the  sidewall.  The  gap  between  the 
rotatable  image  plate  and  the  reflection  plane  shall  be  sealed  with  a sliding  type 
seal  designed  to  minimize  any  flow  leakage.  The  reflection  plane  has  been  test 
calibrated  by  AEDC  with  a movable  nose  flap,  and  acts  as  a boundary  layer  diverter 
which  is  designed  to  channel  the  low  energy  boundary  layer  air  between  the  reflection 
plane  and  the  tunnel  sidewall  while  causing  only  a minimal  distortion  to  the  flow 
pattern  over  the  fuselage. 
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FUSELAGE 


2.  2 

The  fuselage  assembly,  shown  in  GDC  drawing  WT71-107010  #1,  is  the  basic  F-lll  TACT 
1/6'th  scale  fuselage  and  must  be  modified  for  the  forthcoming  test  as  depicted  in  GDC 
wind  tunnel  model  drawing  VVT  7G-1111-08.  Revisions  to  the  fuselage  include 

(i)  provisions  for  internally  mounting  twelve  hydraulic  Moog  valves, 

(ii)  provisions  for  internally  mounting  four  individual  and  one 
(or  more)  manifolded  blocking  valves, 

(iii)  a modified  trailing  edge  glove  fairing  which  permits  approximately  5° 
trialing  edge  deflection  without  fouling  at  the  58°  sweep  condition, 

(iv)  an  extended  leading  edge  fairing  for  testing  with  58°  wing  sweep, 
provided  the  existing  fairing  is  unavailable  or  unsuitable. 

The  model  shall  be  tested  tail  on  and  with  tail  incidence  angles  chosen  so  as 

not  to  exceed  balance  pitching  moment  restrictions.  Tail  settings  shall  be  made  with 

existing  tail  brackets  which  give  incidence  angles  of  0,  -4,  -10,  and  -25  degrees. 

2.  .3  WING  LEADING  EDGE  SECTION 

The  wing  leading  edge  design  consists  of  a front  spar  and  six  independently  controlled 
double  barreled  hydraulic  actuator  systems  as  described  in  GDC  drawing  \VT  7G-111102, 
sheets  1-8.  The  actuator  designations  in  the  leading  edge  section  are  A , A , A 

11  1^  lu 

(inboard)  and  A » A and  A (outboard),  which  vary  respectively  nose  radius, 

^1  wZ  Zu 

leading  edge  deflection  about  the  15‘’e  chord  line,  and  leading  edge  deflection  about  the 
25r,'  chord  line.  The  volumes,  cross-sectional  ai’eas  and  longitudinal  displacements  of 
the  hydraulic  cylinders,  and  the  maximum  angular  and  longitudinal  deflections  of  the  wing 
elements  are  given  in  Table  2. 1.  The  hydraulic  cylinders  for  the  A^?,  A^  and  A^0, 
and  A systems  are  located  inside  the  front  spar  section  between  0.  25c  and  0.40c, 
whereas  the  cylinders  for  the  A and  A systems  are  located  inside  the  nose  section 
(forward  of  0. 15c). 

Potentiometers  for  measuring  and  setting  actuator  position  for  the  leading  edge 

system  are  to  be  mounted  as  shown  in  Figures  2.  2 - 2.4.  Backup  (spare)  potentiometers 

are  provided  for  each  of  the  leading  edge  systems.  For  the  A , and  A systems, 

13  23 
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the  potentiometers  are  linked  between  the  piston  and  the  front  spar  and  measure 
piston  displacement  directly.  For  the  A and  kt>l>  systems,  the  potentiometers  are 
linked  across  the  hinge  line  at  the  15%  chord  to  measure  the  rotation  about  this  position. 
For  the  A and  A ^ actuators,  the  potentiometers  are  mounted  between  the  actuator 
cylinder  wall  and  the  piston,  and  are  designed  to  directly  measure  the  nose  extension, 
provided  the  nose  skin  always  remains  in  contact  with  the  push  rod  which  is  attached 
to  the  pistons.  Experience  with  the  kinematic  model  has  shown  that  this  may  not  be 
the  case  when  the  leading  edge  is  deflected  about  the  15%  chord  hinge  line,  because 
of  the  sliding  skin  joint  on  the  upper  surface  (see  Section  4).  Thus,  deflecting  down 
the  leading  edge  could  unport  the  nose  skin  from  the  push  rod,  due  to  the  lengthening 
of  the  skin.  This  behavior,  if  uncorrected,  would  result  in  large  nonlinear  interactions 
between  the  A and  A^  and  between  A ^ and  Ao  f actuators,  including  considerable 
backlash,  and  could  complicate  the  optimization  procedure  for  nose  radius  by  requiring 


corrections  to  the  calibrations  of  the  A and  A actuators  as  shown  in  Ref.  4, 

11  1 & 

(Eqs.  3. 1 and  3.  2). 

The  leading  edge  skin  is  designed  with  a sliding  joint  (rearward  facing  lap)  on 


the  upper  surface.  The  lap  moves  forward  with  increased  deflection  of  the  A and 

A , and  A and  A actuators.  Care  must  be  exercised  when  mounting  the  leading 
13  22  23 

edge  skin  so  that  the  lap  does  not  import  from  the  underlying  surface  so  as  to  form  a 


gap  which  could  affect  the  data . 


2.  4 WING  TRAILING  EDGE  SECTION 


The  wing  trailing  edge  design  consists  of  the  main  spar  with  four  double  barreled  and 

two  single  barreled  hydraulic  actuator  systems,  as  described  in  GDC  drawings 

WT  76-111103,  sheets  1-4.  The  actuator  designations  are  A,.,  A„„  and  A„„ 

14  15  16 

(inboard),  and  A , A and  A (outboard),  which  vary  respectively  upper  surface 
24  25  26 

humping,  trailing  edge  deflection  about  the  65%  chord  line,  and  trailing  edge  deflection 
about  the  80%  chord  line.  The  volumes,  cross-sectional  areas,  longitudinal  displace- 
ments of  the  hydraulic  cylinders,  and  the  maximum  angular  and  longitudinal  deflections 
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of  the  wing  elements  are  listed  in  Table  2.2. 

The  hydraulic  cylinders  for  the  A.  and  A upper  surface  humping  cylinders 

14  24 

are  circular  in  section  and  are  designed  to  vary  the  spar  thickness  at  the  60%  chord 

position  at  two  spanwise  stations  (SS  36.972  inches  and  54.  167  inches,  corresponding 

to  62.3%  b/2  and  91.3%b/2,  respectively).  The  wing  was  designed  to  simulate  the 

W spar  depth  at  SS  36.972  with  the  A actuator  nearly  retracted,  and  the  W spar 
52  14 

depth  at  54. 167  inches  with  the  actuator  A extended.  The  upper  rear  wing  panel  has 

24 

been  stiffened  with  a spanwise  stringer  along  the  60%  chord  line  to  extend  the  humping 
over  a wider  region  of  the  span.  Additional  descriptions  of  the  various  airfoils  to  be 
approximated  are  given  in  Section  6. 

The  trailing  skins  are  fabricated  in  two  panels  with  sliding  joints  at  the  trailing 

edge,  as  shown  in  Figure  2.5.  The  trailing  edge  thickness  has  been  increased  by  only 

0.25%  c from  that  of  TACT  wing  YV  to  1.  00% c,  necessitating  the  use  of  fasteners  which 

52 

protrude  outside  the  skin  surface  in  the  trailing  edge  region  (Figure  2.5). 

Potentiometers  for  measuring  and  setting  actuator  position  for  the  trailing 

edge  systems  are  shown  in  Figure  2.6  & 2.7.  Backup  potentiometers  are  provided 

for  the  A , A , A and  A . actuators,  but  not  for  the  A and  A humping 
15  16  25  26  14 

systems,  because  of  space  limitations.  In  the  event  of  a potentiometer  failure  for 

either  of  these  two  systems,  the  blocking  valves  should  lock  the  actuator  at  its  last 

position.  The  potentiometers  for  A and  A<?^  read  the  direct  linear  displacement  ol 

the  pistons,  and  hence  the  change  in  spar  box  depth.  The  potentiometers  for  the  A 

and  A actuators  are  mounted  off  of  the  rear  spar,  and  read  the  piston  motion  which 
25 

is  a direct  measure  of  the  rotation  of  the  ti’ailing  edge  about  the  65% c hinge  line.  The 

potentiometers  for  the  A and  A actuators  are  linked  across  the  hinge  line  at  the 

16  26 

80%  position  and  measure  the  trailing  edge  rotation  about  this  hinge  line. 
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TABLE  2. 1 LEADING  EDGE  SYSTEMS 
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TABLE  2.2  TRAILING  EDGE  SYSTEMS 
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Figure  2. 1 Model  Installation  in  Tunnel  16T 
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3.0  SUMMARY  OF  STRESS,  LOADS 
AND  FLUTTER  ANALYSES 


Detailed  stress,  loads  and  flutter  analyses  carried  out  for  the  SOFT  Wing  design  are 
given  in  References  5 and  0.  Some  of  the  basic  assumptions  and  a summary  of  the 
results  are  presented  below. 

3. 1  STRESS  AND  LOADS  ANALYSES 

The  stress  and  loads  analyses  given  in  Reference  5 are  based  on  a free  stream  dynamic 
pressure  qa,  ~ 500  PSF,  and  on  wing  pressure  distribution  data  obtained  from  previous 
wind  tunnel  tests  with  a 1/6'th  scale  semispan  F-lll  TACT  jet  flap  model  (Reference  14). 
The  loads  taken  for  the  stress  analyses  used  the  maximum  differential  pressure 
coefficient  Cpt  - CPu  between  the  lower  and  upper  wing  surfaces  measured  at  each 
tap  on  the  wing.  The  data  were  compiled  from  six  different  runs  covering  wing  leading 
edge  sweep  angles  A.  = 26°  and  45°,  free  stream  Mach  numbers  M*  = 0.  85,  0.  90  and 
0.95.  The  jet  flap  blowing  coefficient  and  jet  efflux  angle  were  Cw=  0.024  and  0j  = 80°, 
respectively.  The  angle  of  attack  a for  the  various  runs  was  selected  as  that  which 
gave  the  maximum  unstalled  lift  coefficient.  Values  of  Cp  ranged  from  0.G5  to  1.48. 

The  resulting  maximum  values  of  CpL  - Cp^  are  plotted  in  Figure  3. 1,  and  were 
used  to  establish  the  wing  pressure  distributions,  forces  and  moments  for  the  stress 
analyses. 

The  design  criteria  called  for  a factor  of  safety  S.  F.  =4  based  on  ultimate, 
and  S.  F.  =3  based  on  yield.  As  seen  from  the  S.  F.  summaries  presented  in  Table 

3. 1  for  the  various  model  details,  these  criteria  were  met  for  the  split  spar  design. 

The  minimum  S.  F.  for  the  spar  is  5.  71  (ultimate)  and  occurs  4 inches  outboard  of 
wing  semispan  station  SS  14.  7.  On  the  other  hand  the  summaries  show  that  various 
actuator  lever  arms  and  brackets  approach  the  S.  F.  =4  limit  very  closely,  with  S.  F . 
values  ranging  from  3.  9 to  4.  G in  ultimate  being  common.  The  page  designations  for 
the  item  or  part  listed  in  Table  3. 1 are  for  Reference  5. 
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Because  of  added  flexibility  requirements,  the  skin  panels  were  designed  solely 
on  the  basis  of  yield,  with  a S.  F.  = 2 being  considered  acceptable.  The  S.  F.  's  given  in 
Table  3. 1 for  the  lower  and  upper  skin  panels.  1.9  and  2.  0,  respectively,  should  be 
considered  approximate  and  are  based  on  zero  differential  twist  between  inboard  and 
outboard  panels.  Similar  values  have  been  predicted  for  the  leading  edge  skin  panels, 
but  in  this  case  the  computations  are  even  more  unreliable,  because  of  the  combined 
effects  of  nose  sharpness  changes,  leading  edge  deflection,  and  airload  on  the  skin 
deformation,  even  without  differential  twist. 

Because  of  the  uncertainty  of  the  stress  levels  in  the  leading  edge  and  trailing 
edge  skin  panels,  it  is  planned  to  implant  approximately  50  strain  gages  in  high  stress 
locations  on  the  outer  surface  of  the  skins  during  the  pretest  calibrations.  These 
stress  levels  will  be  monitored  during  the  calibrations  of  the  actuators,  and  used  to 
establish  limits  for  actuator  deflection  (and  differential  twist)  as  described  in 
Reference  4.  Since  the  gages  must  be  removed  for  the  wind  tunnel  test,  some 
uncertainty  in  maximum  skin  stress  levels  under  airload  conditions  will  remain. 

The  test  operation,  as  discussed  in  Reference  4,  calls  for  making  preliminary  runs 
at  gradually  increasing  values  of  qK,  for  reasons  of  verifying  aerodynamic  damping. 

It  is  planned  to  examine  the  skin  panels  after  each  run  at  increased  q u , in  order  to 
look  for  evidence  of  yield. 

3.2  FLUTTER  ANALYSES 

The  flutter  analysis  of  the  SO I-T  Wing  model,  details  of  which  are  presented  in 
Reference  G,  was  based  on  a highly  simplified  structural  representation  consisting 
of  seven  spanwise  segments.  The  mass,  geometry,  torsional  and  bending  stiffnesses 
for  this  representation,  as  well  as  coupled  bending  and  torsional  modal  and  frequency 
data  are  presented  in  Reference  G. 

The  aerodynamic  damping  computations  were  performed  only  for  the  A = 26° 
sweep  condition  (judged  most  critical).  The  computer  program  was  based  on 
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Theodorsen's  incompressible  unsteady  aerodynamic  strip  theory  model  and  used  a 
lift  curve  slope  value  a = 2 it . The  flutter  equations  were  solved  for  flutter  frequency 
and  damping  coefficient  as  a function  of  air  density  at  a fixed  sound  speed  c = 11G0  ft/sec. 
The  resulting  incompressible  flutter  boundaries  (zero  damping  coefficient)  are  plotted 
in  Figure  3.  2 versus  q^  . Also  shown  in  Figure  3.  2 are  the  compressible  flutter 
boundaries  (obtained  with  the  Prandtl-Glauert  factor)  and  the  envelope  of  proposed 
test  conditions. 

Even  though  the  test  envelope  is  well  within  the  compressible  flutter  boundary, 
various  assumptions  and  uncertainties  in  the  analysis  dictate  that  the  initial  wind 
tunnel  rims  be  made  at  q*  = 300  PSF,  with  qx  increasing  in  steps  of  50  PSF  to  the 
design  value  of  qx  = 500  PSF.  It  is  planned  to  evaluate  the  level  of  aerodynamic 
damping  coefficient  at  each  value  of  q » by  taking  dynamic  data  traces  of  the  response 
of  the  model  (wing  tip  acceleration)  to  step  changes  in  actuator  position,  as  discussed 
in  Reference  4.  It  is  also  planned  to  verify  the  bending  and  torsional  stiffness  values 
used  in  the  flutter  computations  by  making  static  deflection  measurements  during 
model  calibration  (see  Section  5.  2. 1 of  Reference  4). 

Computations  were  also  made  of  the  minimum  skin  thickness  required  to  avoid 
panel  flutter  by  the  semi -empirical  method  given  in  Reference  15.  The  computations 
were  made  for  the  upper  surface,  rear  panel,  which  was  judged  to  be  most  critical, 
and  showed  that  support  was  required  at  the  80$  chord  location.  Because  of  the 
approximate  nature  of  the  panel  flutter  computations,  the  testing  will  be  carried  out 
with  two  strain  gages  on  the  upper  skin  panel  from  which  RMS  readings  will  be 
obtained  on  line  during  the  test. 


TABLE  3.1  SOFT  WING  SAFETY  FACTOR  SUMMARY  FROM  REF. 


TABLE  3.1  (Continued) 


rABLE3.1  (Continued) 


TABLE  3. 1 (Continued) 


TABLE  3.1  (Concluded) 
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Figure  3. 1 Envelopes  of  Maximum  Pressure  Coefficients 
for  use  in  Stress  Analyses 


1.0  EVALUATION  OF  KINEMATIC  MODEL 


A simplified  kinematic  model  was  designed  and  fabricated  to  illustrate  the  operation 

of  the  various  actuation  systems.  The  kinematic  model  was  built  in  two  separate 

parts,  (i)  leading  edge,  and  (ii)  trailing  edge,  as  pictured  in  Figures  4.1-4.  3.  The 

kinematic  models  were  built  to  the  same  scale  as  the  wind  tunnel  model,  but  extended 

only  over  approximately  60',  of  the  wing  semispan,  including  the  semispan  stations 

at  34.  8%  b/2  and  91.  3 7 b/2  at  which  the  airfoil  sections  are  to  be  measured  on  the 

wind  tunnel  model.  The  kinematic  model  of  the  nose  section  was  designed  with  hand 

operated  actuators  for  varying  leading  edge  radius  (A^  and  A actuators),  leading 

edge  deflection  about  the  15','  chord  line  (A  and  A actuators)  and  leading  edge 

12  22 

deflection  about  the  25%  chord  line  (A  and  A actuators).  The  trailing  edge  kinematic 
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model  was  designed  with  hand  operated  actuators  for  deflecting  the  trailing  edge  about 

the  65 % chord  line  (A,  and  A ) actuators  and  about  the  80%  chord  position  (A,  „ and 
la  25  16 

A p actuators).  The  upper  sui’face  humping  mechanisms  (A^  and  A^  actuators)  were 
omitted  from  the  design  to  save  on  cost,  and  because  their  operation  should  be  relatively 
straight  forward  and  not  require  a physical  demonstration.  Separate  actuators  were 
designed  at  the  inboard  and  outboard  stations  of  the  kinematic  model  to  permit 
differential  deflection  across  the  span  of  the  kinematic  model  for  each  of  the  five 
actuator  systems.  This  should  illustrate  the  skin  deformation  and  maximum  stress 
areas  due  to  differential  actuator  deflection  over  that  portion  of  the  skin  simulated  by 
the  kinematic  model. 

The  skin  panels  of  the  kinematic  model  were  fashioned  from  plexiglass  sheet 
stock  which  were  laminated  together  to  approximate  the  thickness  distribution  of  the 
steel  skins  on  the  wind  tunnel  model.  Slip  joints  were  placed  on  the  upper  surface  of  the 
leading  edge  section  and  at  the  trailing  edge.  Thus,  the  skin  deflections  and  stress 
distributions  would  be  similar  to  those  produced  on  the  wind  tunnel  model  with  similar 
actuator  deflections.  In  principal,  at  least,  the  absolute  actuator  force  and  stress 
levels  expected  on  the  wind  tunnel  model  may  also  be  estimated  from  the  kinematic 
model  through  El  scaling  (see  References  10  and  16).  Lack  of  time  and  difficulties 


in  measuring  skin  radius  of  curvature  to  sufficient  accuracy  prevented  this  aspect  of 
the  program  from  being  accomplished. 

An  evaluation  of  the  functionality  of  the  kinematic  model  has  been  given 
in  Reference  1C.  The  kinematic  model  was  also  transported  to  Wright  Patterson 
AFB,  Ohio,  and  demonstrated  at  the  SOFT  Wing  Prefabrication  Meeting  held  there 
on  14  October  197G.  The  kinematic  model  was  also  displayed  and  exercised  at  the 
SOFT  Wing  Pretest  Conference  held  at  GDC  on  13  December  1970.  A summary  of 
the  evaluations  is  given  below. 

4.  1 LEADING  EDGE  UNIT 

Operation  of  the  leading  edge  unit  has  shown  that  the  inboard  and  outboard  actuators 
may  be  deflected  through  the  positions  given  in  Table  4. 1 without  wrinkling  or  deform- 
ing the  nose  skin,  provided  the  same  angular  deflections  are  maintained  inboard  and  out- 
board, i.e.,  zero  differential  twist  or  camber.  However,  considerable  crinkling  and 
warping  of  the  nose  shape  was  observed  whenever  the  total  inboard  deflection  angle 
(as  given  by  A + A ) differed  from  the  total  outboard  deflection  angle  (as  given 

by  A + A ) by  as  much  as  10°. 
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TABLE  4.1  LEADING  EDGE  DEFLECTIONS 


INBOARD 

OUTBOARD 

A13  (Deg) 

A12  (Deg) 

A23  (DCg) 

A22  (Deg> 

20 

0 

20 

0 

20 

10 

20 

10 

20 

20 

20 

20 

10 

0 

10 

0 

10 

10 

10 

10 

10 

20 

10 

20 

0 

0 

0 

0 

0 

10 

0 

10 

0 

20 

0 

20 

J 
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A second  potential  difficulty  encountered  in  the  nose  region  stemmed  from  the 


bending  of  the  rod  which  is  pushed  into  the  nose  of  the  airfoil  by  the  A and  A nose 
sharpness  actuators.  Although  the  rod  is  much  stiffer  in  the  kinematic  model  (made 
of  aluminum)  in  relation  to  the  skin  than  in  the  wind  tunnel  model,  considerable 
bending  of  the  rod  in  the  chordwise  direction  w'as  observed  where  it  is  unsupported 
by  the  actuators.  The  nose  sharpness  push  rod  in  the  wind  tunnel  model  will  be 
supported  in  additional  locations  and  stiffened  considerably  to  prevent  chordwise  bend- 
ing from  occurring. 


A third  difficulty,  which  was  observed,  dealt  with  the  unporting  of  the  nose 
sharpness  push  rod  from  the  nose  skin  when  the  A and  A actuators  are  displaced 
to  deflect  the  nose  about  the  151  chord  position  while  holding  the  A and  A nose 
sharpness  actuators  fixed.  If  uncorrected,  this  would  make  the  variation  of  nose 
radius  highly  nonlinear  with  regard  to  the  displacement  of  the  A and  A actuators, 
lead  to  a sizable  back  lash  zone  over  which  leading  edge  radius  does  not  vary  with 


actuator  position,  and  lead  to  interactions  between  the  A^  and  Ao^  actuators  and  the 

A and  A systems.  Aerodynamic  loads  on  the  skin  might  pull  the  nose  skin  up 
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tightly  around  the  push  rod  under  the  conditions  of  high  lift  and  high  tunnel  dynamic 


pressure,  making  the  actual  variation  of  nose  radius  dependent  on  these  tunnel 


parameters  as  well  as  upon  the  A , A , A , and  A actuators. 

XX  uX  X^  uu 

One  possible  method  of  avoiding  the  difficulty  may  be  to  program  the  A^  and 

A actuators  to  displace  linearly  with  the  position  of  the  A and  A systems.  This 
21  XZ 

would  require  additional  pi'ogramming  changes  to  both  the  optimization  program  and  to 
the  controller,  as  shown  in  Reference  4.  However,  a final  decision  in  regard  to  this 
interaction  will  not  be  made  until  the  wind  tunnel  model  has  been  calibrated  and  the 


same  effects  observed. 


4.2  TRAILING  EDGE  UNIT 


Initial  operation  of  the  trailing  edge  unit  has  shown  that  the  inboard  and  outboard  actuators 
may  be  deflected  through  the  positions  shown  in  Table  4.  2 without  w rinkling  either 
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the  upper  or  lower  skin.  In  contrast  to  the  observations  made  with  regard  to  the 
leading  edge,  no  crinkling  or  undue  warpage  of  the  skin  was  produced  even  with  a 
10  difference  in  deflection  angle  between  the  inboard  and  outboard  systems, 
although  a bulge  was  noted  at  the  trailing  edge  for  this  condition,  and  excessive 
stress  may  be  introduced  in  the  skin  panels  aei'oss  the  spanwise  gap  between  the 
inboard  and  outboard  actuators.  It  is  pi  aimed  to  implant  strain  gages  on  the  skin 
of  the  wind  tunnel  model  in  the  gap  region  to  check  the  stress  levels  produced,  and 
to  increase  the  span  of  the  gap  between  the  actuators,  should  this  prove  necessary. 

TABLE  4.2  TRAILING  EDGE  DEFLECTIONS 


INBOARD 


— 

A 

15 

A1C 

20 

0 

20 

10 

20 

20 

10 

0 

10 

10 

10 

20 

0 

0 

0 

10 

0 

20 

OUTBOARD 

A 

25 

2f» 

20 

0 

20 

10 

20 

20 

10 

0 

10 

10 

10 

20 

0 

0 

0 

10 

0 

20 

5.0  COMPUTER  SIMULATION  STUDIES 


The  optimization  computer  program  plays  a central  role  in  the  SOFT  Wing  wind  tunnel 
testing  technique  being  planned  at  AEDC,  in  that  the  optimization  program  determines 
how  to  change  the  actuator  settings  and  angle  of  attack  in  order  to  achieve  the  optimum 
value  of  a specified  aerodynamic  function  subject  to  given  constraints.  The 
relation  of  the  optimization  program  to  the  SOFT  Wing  wind  tunnel  testing  operation 
is  shown  in  the  flow  chart  in  Figure  5.1. 

The  optimization  program  has  controlling  parameters  which  specify  which  aero- 
dynamic parameter  is  to  be  the  objective  function  and  which  of  the  remaining  parameters 
are  to  be  subject  to  equality  or  inequality  type  constraints.  The  output  from  the  optimiza- 
tion program  is  an  array  of  actuator  settings  and  angle  of  attack.  The  wing  is  subsequentl 
positioned  according  to  these  settings.  Wind  tunnel  data  are  then  taken  and  the 

results  reduced  to  coefficient  form,  e.  g.  C , C , C , etc.  These  reduced  wind 

L D M 

tunnel  data  are  then  input  to  the  optimization  program,  so  that  it  can  evaluate  the 
effect  of  its  previously  predicted  settings.  New  settings  are  then  calculated  and 
output  in  an  attempt  to  optimize  the  objective  function. 

The  SOFT  Wing  testing  technique  is  basically  a nonlinear  programming 
optimization  procedure  in  which  the  optimization  computer  program  views  the 
wind  tunnel  as  an  analog  computer  for  evaluating  the  following  mathematical 
functions:  CL(Aij,rv) 

C d = Cp)  (A. Oi  ) 
ij 


CM  Cm  (A. . , oc ) 


'M 


ij 


(3.1) 


c 


P-r 


TE1,  2 
etc. , 


= C 


P-, 


TE 


1,2 


<V> 
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where  these  and  other  possible  aei'odynamic  parameters  are  defined  in  the  Nomenclature 
Section.  Thus  the  wind  tunnel  system  consisting  of  the  airstream,  test  section,  support 
system,  model,  control  system,  instrumentation  and  data  reduction  package,  is  considered 
simply  as  a physical  function  generator  with  which  the  optimization  program  interacts. 

The  advantages  of  using  the  wind  tunnel  as  a function  generator  rather  than  making 
analytical  computations  of  the  aerodynamic  parameters  during  the  optimization,  as 
done  in  References  17  - 19,  stems  from  the  improved  accuracy  and  speed  of  the  wind 
tunnel  over  any  known  flow  field  computational  procedure.  This  will  be  especially  true 
under  conditions  of  viscous,  turbulent,  transonic  flow  and  for  complicated  three 
dimensional  geometries  representative  of  aircraft.  The  disadvantages  of  using  the 
wind  tunnel  stem  from  wall  type  interactions  and  scale  effects,  which  may  make  the 
results  difficult  to  apply  to  the  flight  of  full  scale  aircraft,  to  the  restrictions  imposed 
by  limits  of  model  articulation  and  to  instrumentation  errors. 

Section  5. 1 presents  the  mathematical  basis  for  the  optimization  computer 
program  to  be  used  during  the  forthcoming  test.  Section  5. 2 presents  a highly 
simplified  set  of  algebraic  equations  used  to  calculate  the  aerodynamic  coefficients 
during  computer  simulation  studies  of  the  optimization  procedure.  These  algebraic 
equations,  termed  "Pseudo  Aerodynamic  Module,"  replace  the  entire  wind  tunnel 
system  when  conducting  test  simulations  with  the  optimization  program.  Section 
5.  3 presents  a number  of  sample  results  from  the  simulations  illustrating  the  effects 
of  various  parameters  in  the  optimization  program  on  convergence  rate.  As  a result 
of  the  simulation  studies  it  was  concluded  that  all  optimization  problems  to  be 
attempted  during  the  wind  tunnel  test  should  be  first  simulated  by  using  the  Pseudo 
Aerodynamic  Module,  in  order  to  establish  the  proper  values  for  the  various 
optimization  parameters  prior  to  the  test. 


5.1  OPTIMIZATION  PROGRAM 

The  present  section  describes  the  mathematical  basis  for  the  optimization  computer 
program,  and  presents  the  input  parameters  which  define  the  wing  optimization 


problem  to  be  solved  and  control  the  rate  of  convergence.  Additional  programming 
details  are  given  in  Reference  4. 

5. 1. 1 Mathematical  Basis 

The  optimization  program  uses  the  gradient  projection  method  for  numerically  solving 

the  constrained  optimization  problems  arising  from  wing  optimization,  hi  conjunction 

with  this  procedure,  transformations  of  variables  are  used  to  enforce  independent 

variable  constraints,  and  restoration  steps  are  required  to  restore  equality  constraints 

that  have  become  unsatisfied  due  to  nonlinearities. 

The  actuator  readings  A..,  tail  incidence  angle  oc  and  wing  angle  of  attack 

i]  t 

a , which  constitute  the  independent  variables  in  the  physical  coordinates,  are 
represented  internally  in  the  program  as  the  vector  0 , where 

A (XMAX)  sin“  ( 0 ) (5.2) 

ij  n 

a 2 (AT  MAX)  sin“  9 ^ - AT  MAX  (5.3) 

2 

and  a 2 (AMAX)  sin“  0 - AMAX  (5.4) 

14 

Here  the  subscripts  i,  j and  n have  the  correspondanee 
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whereas  XMAX,  ATMAX  and  AMAX  are  parameters  which  specify  the  maximum  actuator 
readings,  maximum  (and  minimum)  tail  incidence  angles,  and  maximum  (and  minimum) 
wing  angles  of  attack. 

The  optimization  procedure  used  by  the  optimization  computer  program  is  bas  > 
upon  the  gradients  of  the  objective  function  and  the  equality  constraint  functions.  The 
independent  variables  are  changed  so  that  the  objective  function  decreases  in  a modified 
steepest  descent  direction.  The  direction  is  modified  by  projecting  the  negative  gradient 
of  the  objective  function  onto  the  tangent  planes  of  the  equality  constraints.  In  this  way 
the  independent  variables  are  changed  so  that  the  design  moves  along  a linear  projection 
of  the  constraints.  If  g is  a row  vector  of  the  equality  constraint  functions,  all  of  whose 
components  are  zero  when  the  equality  constraints  are  satisfied,  0 is  the  objective 
function  to  be  minimized,  F is  the  Lagrangian  function,  and 

!+l  if  0 is  to  be  minimized 

(5.5) 

-1  if  0 is  to  be  maximized 
then  the  search  direction  used  is 


s = -(C  V 0 + V g A ) = -V  F 

(5.6) 

where 

A = -C  ( V gT  V g)-1  VgTV0 

(5.7) 

and 

F = C 0 +gA 

(5.8) 

The  independent  variables 

are  changed  by  adding 

a 

ii 

< 

(5.9) 

where  K is  the  step  magnitude  which  is  selected  to  optimize  the  Lagrangian 
function  (5.  8)  in  the  direction  S. 

Although  the  gradient  projection  method  gives  the  search  direction  S,  it  does 
not  yield  the  magnitude  of  the  step  size  to  be  taken.  This  must  be  found  by  taking  a 
series  of  steps  in  the  search  direction  and  then  determining  which  step  size 
minimizes  F.  Since  all  independent  variables  are  incremented  simultaneously  during 
this  process,  this  operation  has  been  termed  the  "simultaneous  mode"  portion  of  an 
optimization  run.  That  portion  of  the  optimization  procedure  during  which  the 
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independent  variables  are  perturbed  individually  and  then  reset  to  their  basic  position, 
in  order  to  generate  the  components  of  the  gradient  matrices  a 0 and  7 g , has  been 
termed  the  ''incremental  mode"  portion  of  the  optimization  run  (See  also  Hefs.  4,  9, 

10  and  11). 

Eleven  values  of  K are  used  during  the  simultaneous  mode.  They  are  selected 
as  follows: 

Ky  = 0 

K2  = the  smallest  value  such  that  one  of  the  components  of 
^ 0 is  two  times  the  standard  deviation  of  error  in 
that  component. 

Kk  = Kk-1  + (1.  4)k~“  K.-,  k = 3,  4 11 

If  K*  denotes  the  particular  Kk  at  which  the  Lagrangian  is  minimal  then  the 
new  independent  variables  for  the  next  incremental  mode  are 


new 


old 


K*  S 


old 


(5.  10) 


and  the  process  is  repeated  if  no  restoration  is  required  The  accomplishment  of 
Equation  5. 10  during  the  wind  tunnel  test  requires  that  all  actuators  and  wing  angle 
of  attack  (tail  incidence  angle  will  be  held  fixed),  be  reset  with  a high  degree  of 
precision  to  avoid  compromising  convergence. 

The  process  of  finding  K*  during  the  simultaneous  mode  operation  is  called 
a one  dimensional  search.  For  each  new'  step  Kk  oi  the  one  dimensional  search,  all 
constraint  functions  are  checked.  If  any  equality  constraint  (or  activated  inequality 
constraint)  drifts  too  far  from  being  satisfied,  then  the  one  dimensional  search  is 
aborted.  If  any  inequality  constraint  becomes  violated  beyond  a prescribed  limit 
during  the  one  dimensional  search,  the  search  is  also  aborted.  In  the  event  of  an 
abort,  the  point  which  had  the  lowest  value  for  equation  (5.  8)  but  did  not  cause  the 
abort  is  taken  as  the  optimal  point  for  the  one  dimensional  search. 

The  components  of  the  vector  DEV  (see  Section  5.1.2)  specify  the  absolute 
values  of  the  maximum  deviations  permitted  for  all  constraints  without  aborting  the 
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the  one  dimension  search.  The  components  of  the  vector  TOL  (Section  5.1.2)  specify 
the  tolerance  to  which  the  constraints  must  be  satisfied  prior  to  searching  in  the 
optimization  direction  S. 

Restoration  is  the  process  of  finding  a point  0 at  which  all  equality  and  activated 
inequality  constraints  are  satisfied  within  TOL.  As  noted  above,  the  one  dimensional 
search  can  result  in  a point  0 at  which  one  or  more  constraints  have  drifted  beyond 
TOL,  due  to  the  nonlinear  nature  of  the  aerodynamic  functions  (dependent  variables). 
When  this  occurs,  a restoration  run  must  be  carried  out.  The  search  direction  which 
will  accomplish  this  (in  the  linear  sense)  is 

S =Vg  (Vgr  V g )-1  g (5.11) 

r 

As  in  the  optimization  case,  a one  dimension  search  must  be  made  to  determine  the 
K*  such  that  the  restoration  criterion  function 

R-RgT  (5.12) 

is  minimized.  However,  in  the  case  of  restoration  with  linear  dependent  variables 
the  theoretical  stepsize  for  minimizing  R is  K*  =1.  Therefore  the  candidate 
stepsizes  are  selected  as  follows  for  restoration: 
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Figure  5.  2 summarizes  the  flow  of  the  mathematical  operations  just  described. 
Each  time  a new  point  0 is  determined,  all  inequality  constraints  are  checked  to 
determine  whether  they  are  still  satisfied.  If  some  are  not,  then  those  that  arc  not 
arc  activated,  which  means  they  are  added  to  the  set  ol  equality  constraints  which  must 
be  kept  satisfied  by  the  operations  described  in  the  proceeding  section.  Similarly, 
those  equality  constraints  which  are  the  result  of  an  activated  inequality  constraint 
on  a previous  step  must  be  checked  to  determine  whether  they  can  be  removed  from 
the  set  of  equality  constraints.  The  condition  for  removal  is  that  the  current  search 
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direction  must  be  interior  to  the  region  defined  by  the  inequality.  Mathematically,  this 
condition  for  constraint  & is, 


i 

j 


(VF)T  VgjL  <;  o.  (5.  13) 

If  this  condition  holds  the  equality  constraint  g ^ is  deactivated. 

A procedure  which  is  employed  to  improve  the  convergence  rate  ol  the 
gradient  projection  method  is  the  Davidon-Fletcher-Powell  (DFP)  variable  metric 
method.  This  procedure  differs  from  the  gradient  projection  method  previously 
described  by  the  inclusion  ot  a scaling  matrix  II  and  by  formulas  for  updating  11. 

For  this  method,  equations  (5.  6)  and  (5.  7)  become 

S = -H  (C  V0  + Vg  A)  (5.  14) 

A = -c  (Vg^  IlT’g)"1  Vg  1 H V 0 (5.15) 

The  matrix  H is  initially  the  identity  matrix  in  which  case  equations  (5. 14)  and 
and  (5. 15)  reduce  to  equations  (5. 6)  and  (5.7).  When  two  consecutive  minimization 
iterations  occur  (without  an  intervening  restoration)  second  derivative  information 
can  be  computed  by  differencing  first  derivative  values: 

A (VF)  = V F -vf0ih  (5.  If-) 

' ' new  oia 


If 


A 0 = -K*  II  V F 


(5.17) 


then  the  update  formula  for  H is 

+ _A8A91 Hold  A(VF)A(7F)Tuold  (5.18) 

"new  old  A0TA(7F)  A ( V F)T  H 0ld  A (v  F ) 

For  nonquadratic  functions  F,  there  is  no  guarantee  that  ll,lew  will  be  an  improvement 
over  H0ic]>  consequently  1^^  must  be  checked  for  validity . Since  11  must  be 
positive  definite,  it  must  have  positive  diagonal  elements.  A check  for  positive 


f 


j 
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diagonals  is  performed, and  if  any  negative  values  are  found  the  Hnew  is  reset  to  the 
identity.  The  other  conditions  for  which  Hnew  is  reset  to  the  identity  are  when  K*  = Kj 
in  the  one  dimensional  search,  and  whenever  a restoration  is  performed. 

5.1.2  Input  Parameters 

Prior  to  attempting  a particular  vving  optimization  problem,  either  in  the  wind  tunnel 
or  in  a simulation,  the  problem  to  be  solved  must  be  specified.  This  is  accomplished 
by  inputing  to  the  program  an  ordered  list  of  up  to  18  aerodynamic  functions  such  as 
shown  in  Equation  5. 1.  These  functions  are  recognized  by  the  program  as  the  array 
FX(I)  of  dimension  I.  To  specify  which  of  the  18  functions  is  to  be  optimized,  the 
variable  IOBJ  must  be  set  to  an  integer  value  from  1 to  18. 

To  define  equality  and  inequality  constraints  the  array  CAN(I),  1=1,  18  must 
be  set  If  CAN  (I)  - 0,  then  function  I is  not  used  at  all.  If  CAN(1)  = 1,  then  function 
I is  an  equality  constraint.  Finally,  if  CAN (I)  = -1,  then  function  1 is  an  inequality 
constraint. 

The  vector  of  values  T(I),  1 = 1,  18  are  the  target  value  for  the  constraints. 

For  equality  constraints,  FX(I)  = T(I)  is  the  mathematical  condition  constraining  the 
optimum.  For  inequality  constraints  FX(I)  a T(I). 

The  angle  of  attack  a is  bounded  by  the  variable  AM  AX  which  is  nominally 
set  at  20  degrees:  | a AMAX. 

The  actuator  positions  are  bounded  by  the  variable  XMAX  which  is  nominally 
set  at  1000  counts:  0 ^ Ajj  s XMAX. 

The  initial  positions  of  the  actuators  and  the  initial  angles  of  the  tail  and  wing 
are  specified  by  the  vector  AVE  (N),  where  the  correspondance  between  N and  the 
independent  variables  is  as  given  in  the  table  subsequent  to  Equation  5.4. 

During  the  wind  tunnel  test  all  independent  variables  may  not  be  used  (active) . 
Thus  at,  corresponding  to  N = 13,  will  be  fixed  on  the  model  during  an  optimization, 
since  the  model  does  not  have  a self  trimming  tail.  Also,  various  actuators  will  be 
deactivated  to  simulate  the  reduced  articulation  of  the  AFTI-111  aircraft 
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(See  Table  6. 1 of  Reference  4).  Finally  various  actuators  may  fail  during  the  course 
of  the  wind  tunnel  t >st  (blockage  valves  in  the  model  should  permit  testing  with  the 
remaining  actuators,  should  either  the  A^,  A ^ and  A ^ or  A systems  fail).  To 
provide  for  the  deactivation  of  an  independent  variable  the  array  DA(N)  is  used.  If 
DA(N)  0,  then  AVE  (N)  is  the  frozen  value  of  independent  variable  number  N. 

The  vector  DA(N)  is  also  used  to  specify  the  magnitude  of  the  perturbations 
during  the  incremental  mode  operation.  If  DA(N)  ^ 0,  then  independent  variable  N > 
is  perturbed  by  DA(N),  the  direction  of  the  perturbation  being  away  from  the  closest 
boundary  (see  Reference  4). 

ZNL  (I)  is  the  vector  with  components  equal  to  the  one  standard  deviation 
noise  levels  (testing  errors)  for  the  measured  values  of  the  independent  variables  I 
14  and  the  18  function  values  fori  15,  1C,  ...32.  The  ZN  L vector  is  utilized  in 
determining  the  size  of  the  K?  step  during  the  simultaneous  mode  portion  of  an 
optimization  run,  as  discussed  in  connection  with  Equations  5.9  and  5. 10.  The 
ZNL  vector  is  also  used  to  introduce  effects  of  testing  type  errors  into  the  analytical 
simulations,  as  will  lie  shown  in  Section  5.3. 

As  noted  previously  in  Section  5. 1. 1 in  describing  the  one  dimensional 
search  and  restoration  procedures,  TOL(I)  is  the  tolerance  for  each  of  the  18 

aerodynamic  function  values  which  are  equality  or  activated  inequality  constraints.  3 

The  constraint  FX(I)  is  assumed  satisfied,  or  restored,  if 

|FX(I)  - T(I)  j £ TOL(I)  (5.19) 

DEV(I)  is  the  tolerance  used  during  the  one  dimensional  search  to  determine 
whether  an  equality  or  inequality  constraint  becomes  violated.  This  is  taken  larger 
than  TOL(I)  to  avoid  restricting  progress  of  the  one  dimensional  search  while  at  the  same 
time,  not  allowing  the  constraints  to  be  violated  by  too  much.  If  either 

FX(I)  - T(I)  > DEV(I)  (for  equality  constraints  (5.  20) 
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or 


FX(I)  - T(I)  + TOL(l)  • 0 (for  inequality  constraints)  (5.  21) 


then  the  one  dimensional  search  is  aborted. 

CMAX  specifies  the  number  of  times  cycled  through  the  incremental  mode. 
In  the  presence  of  noise  (testing  errors),  redundant  perturbations  of  the  independent 
variables  are  taken  and  the  results  are  averaged  together  to  improve  the  accuracy  of 
the  gradient  components. 

A summary  of  all  input  parameters  is  included  in  the  dictionary  of  common 
variables  listed  in  Table  5. 1. 

5. 1.3  Tabular  Input  for  Simulating  a Self-Triming  Tail 

As  noted  previously,  the  optimizations  to  be  carried  out  with  the  present  model 
must  be  made  with  the  horizontal  tail  held  fixed,  since  a variable  tail  incidence 
mechanism  was  not  incorporated  in  the  model.  Hence,  the  constraint  on  trim 

CM 

can  not  be  implemented  directly. 

An  approximate  method  of  implementing  the  trim  constraint  has  been 
presented  in  R Terence  20.  The  method  introduces  two  inequality  constraints  to 
insure  that  sufficient  tail  trim  control  is  available  for  all  wing  shapes  reached  during 
the  optimization,  viz. , that 


CM  ) S CM  (a , at*.  V) » - CM  Mra «“  • “r  V 


(5.  22) 


M ' 1 ’ XJ'  "M 

and  - CM  (or,  , A..)  * CM  ^ (a , <*,  Aj*  ) -CM  V (5-  23) 

Here  C (a , a! , A*  A is  the  pitching  moment  coefficient  as  measured  during 
M 1 k)  + 

the  wing  optimization  runs  at  a fixed  tail  incidence  ~ Qt  an<J  at  variable  or  and 

* * 

Ajj's.  C]yj  (or,  or  j.  , Ay ) is  the  pitching  moment  variation  with  or  recorded  during  a 
preliminary  run  at  fixed  tail  angle  or^  and  with  the  actuators  held  fixed  at  the  nominal 
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positions  A^  A^j.  CM  (a,  at,  Ay ) represents  the  variation  with  a of  the 
minimum  pitching  moment  coefficient  as  recorded  during  a series  of  preliminary  runs 
at  different  tail  incidence  angles  at  and  at  the  nominal  actuator  settings  Ajj  . Cm  MAX 
(Q  , o- 1,  Ajj  ) represents  the  corresponding  variation  for  the  maximum  pitching  moment 
coefficient. 

The  coefficients  (a,  0^,  Ajj  ) and  Cm  max  ^ ’ 01 1’  ^ mus*'  ,)e 

input  in  tabular  form.  The  tables  are  designated  in  the  optimization  program  as 

ALP(1)  « 

CM  MIN  (I)  = CM  min  ( a - Aij  > 

CMNOM(I)=CM(o,a*t.  Ajj) 

CM  MAX (I)  =cMmAX(“’ AiJ  ) 

LTAB  = number  of  I (points  in  ALP  array). 

It  is  assumed  that  the  tables  are  functions  of  the  independent  variable  array 
ALP.  A separate  set  of  tables  is  generally  required  for  each  Mach  number. 

When  using  the  tabular  form  of  the  trim  constraint,  it  is  meaningful  to 
correct  Cj  for  the  change  in  lift  which  would  have  occurred  had  the  tail  been 
self  trimming.  The  correction,  as  derived  in  Reference  20,  is 


= Hr)-1  CM(a.  - Aii) 
i'RTTU  ' c ' M 1 iJ 


lTRIM 


(5.24) 


and  depends  upon  the  ratio  of  the  tail  arm  to  Ule  mean  aerodynamic  cord  c and 
the  measured  Ct^j.  The  latter  coefficient  depends  upon  the  longitudinal  position  of 

the  center  of  gravity  (or  moment  reference  center  x Mc ) • Moving  the  C.G. 
forward  makes  CM  more  negative  and  detx-acts  from  Cl.  whereas  moving  the  C.G. 
aft  has  the  opposite  effect.  Hence,  the  optimization  procedure  may  be  used  to 
investigate  effects  of  static  margin  on  wing  optimization. 

The  corrected  lift  coefficient 


cL  =cl+acLtrim 

must  then  be  used  as  one  of  the  eighteen  aerodynamic  functions  in  FX(I). 


(5.  25) 


4H 


5.2 


PSEUDO-AERODYNAMIC  MODU L K 


Computer  simulations  are  required  to  check  the  convergence  of  the  optimiza- 
tion procedure  prior  to  testing,  because  of  the  various  nonlinearities  and  testing 
inaccuracies  which  could  seriously  compromise  convei'gence.  This  is  especially  true 
under  transonic  conditions  and  with  a large  number  of  independent  variables,  since  these 
effects  tend  to  increase  the  sensitivity  of  the  computed  optimum  search  directions  S 
to  testing  errors,  and  make  the  restoration  and  the  holding  of  constraints  increasingly 
difficult. 

hi  order  to  carry  out  computer  simulations  of  the  interactive  testing  operation, 
it  was  necessary  to  link  the  optimization  program  with  a set  of  equations  for  calculating 
the  aerodynamic  functions  ordinarily  supplied  by  the  wind  tunnel  data  reduction  system 
(see  Equation  5. 1).  These  equations,  termed  "Pseudo-Aerodynamic  Module"  (PAM) 
are  a set  of  nonlinear,  algebraic  relations  for  the  aerodynamic  functions  of  Equation 
5. 1 in  terms  of  the  twelve  actuator  positions  Ajj,  Aj2>  • • • ^26’  incidence  angle 
0%,  angle  of  attack  a,  Mach  number  M and  Reynolds  number  Re.  Because  of  the 
numerous  times  that  the  PAM  subroutine  must  be  called  upon  during  an  optimization 
problem,  it  was  impractical  from  a computer  time  standpoint  to  perform  even  the 
most  simplified  type  of  flow  field  analysis.  Hence  emphasis  in  PAM  was  placed  on 
computational  speed  and  on  a realistic  approximation  of  the  form  of  the  nonlinear 
features  expected  in  the  aerodynamic  test  data,  rather  than  on  accuracy.  Random 
errors  were  also  added  to  the  ouput  from  PAM  to  simulate  testing  inaccuracies.  The 
free  stream  Mach  number  M in  PAM  must  be  taken  as  subsonic,  viz. , M < 1. 

The  following  expressions  constitute  PAM  as  currently  programmed,  and 
are  based  somewhat  on  thin  airfoil  theory  with  empirical  corrections  for  three 
dimensional  and  compressibility  effects.  Numerical  values  for  the  various  coefficients 
were  chosen  to  correspond  with  the  1/6'th  scale  SOFT  Wing  wind  tunnel  model. 
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5.2.  1 Relation  Between  Actuator  Readings  and  Geometric  Deflections 
(R/c  ) 0.01  (3-2A11/l000  ) (inb'dL.E.  radius) 

6 25  (A^/1000  - 1/5)  (inb'd  L.E.  deflection,  degs,  15r,c) 

6 ^ 25  (Ari/1000  - 1/5)  (inb'd  L.E.,  degs,  25%c) 

(t/c)  0.09  + 0.01  (2A  /1000  - 1)  (inb'd  upper  surface  humping) 

14  14 

6 30  (A,  / 1000  - 1/3)  (inb'd  T.E.  deflection,  degs,  65%c) 

15  15 

6 30  (A^/  1000  - 1/3)  (inb'd  T.E.  deflection,  degs,  80l£c) 

(R/c) § 0.  005(  3 - 2 A _ /1000)  (outb'dL.E.  radius) 

6 ( 25  (A  } /1000  - 1/5)  (outb'd  E.E.  deflection,  degs,  15%c) 

f)  25  (A^/1000  - 1/5)  (outb'd  E.E.  deflection,  degs,  25r<c) 

(t/c),  O.Ofi  +0.01  (2  A /1000  - 1)  (outb'd  upper  surface  humping) 
24  24 

6 =30  (A./1000  - 1/3)  (outb'd  T.E.  deflection,  degs,  65%c) 

25  25 

6 = 30  (A^/lOOO  - 1/3)  (outb'd  T.E.  deflection,  degs,  80%c) 

2G  26 


5.2.2  Aerodynamic  Coefficients 

(i)  Section  Zero  Lift  Angle 

"ol i'TT  A(CL6)i -“TWj 


, i =1,  2 


<0^  = 2 £(sin  02  - 02)sin6. + (sin  03  - 03)  sin 6 . 3 

+ T)i5(sin  05  + 05)  sin  6 _5  + ni6(sin  0fi  + 0f))  sin  6i(?j 


i = 1,  2 


0 j = cos  [l  - 2 Cj/c  ] 
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and 

^15  =T)16  =T)25  =T126  = 1- 

(ii)  Wing  Lift  Curve  Slope 

v Al 

cha  = - —UT 

1 + 1 1 + [l-  (M  cos  jV)2](A1  /2  cosjv)  | 

where  the  values  At  =5.83,  -A=  16°  and  M = 0.  85  were  input  for 
the  simulations. 

(iii)  Section  Stall  Parameters  and  Section  Lift  Coefficient 

i/3  1/2  1/2 

Clmax.  =0-8 +0-3  <R/c)ii  +4(c2/c)  Sin5i2 +4(c3/c)  sin6i3 

+ % (sin  05  + fl5)  sin6i5  +lli6(sine6  +06)sin6i6 
a STALL  j =a0Li  +57'3  CLMAXi/CLa 
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STALL a 





r a 

Ci  . = Cr  1-H(- 
Li  L a 


-1)]  sin  (“-“oLj' 

H (_ S 1)CL  sin  STALL i ' a0Li>  [*  "H(cTr*"1'] 

a STALL  i ' a 1 L ci 


T “ _a  STALL  i .1  „ a _1X 

-AC,  . 1-  cos  90(  — ) H(  i) 

lL  aC,  '“STALL,-1  “STALL. 


H'rr1*] 


a 

1.  8 sin(a-a0L  ) cos  (a-a01  ) H (— 

1 ^ i 


-1) 


where 


aci  “stall. 


+ [lOO  (R/c)41]' 


A CL  = sin 


^“c  j ~ “ OL  [CLa 


sin  (“STALL  ^“OLi*  _ x 

sin(ttc. -aoLj) 


8 COS  (Q! 


and  H(x)  is  the  unit  step  function,  H=0,  x<  0,  H 1,  x > 0. 

(iv)  Component  and  Total  Lift  Coefficient 
a.  Wing 


lWING 


AK1  CL i +^K2  CL2 


with  A K i = 0.41  and  A K 2 - 0.  24 


b.  Glove 


°L GLOVE 


=,  (1  - AKi  - A K 2 ) CL^  sin»  cosa 
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c.  Tail 


= 2. 12  and 


is  given  by  the  same  expression  as  C j except  use 
At  - 57°  instead  of  the  corresponding  wing  values. 


The  downwash  angle  e t (in  degrees)  is 
90  ffLVVEO  +°L  GL0V£  ^yi+(0.785  b/<t)2] 


(0.  785  V)2 


where  the  span  to  tail  ax-m  ratio  h/ ^ - 2. 5. 


Lt=CLatsin(Q;t-et)  f1-^  — 


a + Q?t  - et 


STALLt 


- 1) 


* CLo,t  “‘“'"stall,1  h < aaJALl~ 11 

r /«  +Qt~et  ~gSTALLt\l  la  +Q,t~et  _,\ 

-ACLt  L1"0089^  ®ct  ■ ^STALLt  'J*'  aSTALLt  ' 

[ of  +<**.  - e+  , 1 

x 1 - H ( i i 1) 

L . J 


la+a  t“et  \ 

H 1 

f V ' 

t “ct 


-AC, 


where 


(a  +Qft  - et) 


a STALL  t " 12  a +at  -et  ‘ 


(in  degs) 


(<*  +at  - et> 
ct  a + - e ^ 


=14 


(in  degs) 
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and  A C 
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1 


(vi)  Wing  Compressibility  Drag 


AcD  =PLi(w  (MMAXi  -1-05)  H(MMAXi-1-°5)[l-«(MMAXf1-45)] 


+ (MMAXj  -1-2)2  11(MMAXj  " 145)  \ ’ 


i = 1,  2 


vvhe  re 


10 


p =10 

L9 


10  (t/e)2  +100  (R/c)“  -y  (°^  .j_- 

14  11  4 57.3  I 

16.7(t/c f + 200  (R/c)2 


21 


57.3 


and 


M 


1/2  r 


MAXj  Vy 


y-i  9 y M 
(l  + — y m )(  2 cp 


M IN  j 


1) 


Jlzi) 

i 1 y '-i 


1/2 


Here  CTi  is  the  minimum  pressure  coefficient  reached  on  the  i'th 

p MINj 

section  of  the  wing,  and  is  taken  as  the  more  negative  value  of  the  following  two 
expressions,  but  not  less  than  the  third  expression,  which  is  related  to  the  vacuum 
pressure  coefficient. 

jl_  r 

, s - r. 

’MINi 


'P*!  -yTTS  j2(t/Oi4+10(R/c)ll+0.»5CIjt(rt«(«-<»0L1) 


or 


c ' °-2  [cLtt  aln(tt-otoLi>j  _ sin0  sinfi.y  sin 02 sin6i2 J 

ImINj  / 1-M2  I (t/c)-4  + 5 (R/c)jj 


but 


C 


1.5 


pMINi-  yM2 


(vii)  Total  Drag  Coefficient 


CD=CD  + C 


Dtv 


F ™ WING 


aki(ACd„  4ACdrl.) 


1 


AK2(ACd  +acd 
a o 


RL . 


+ ^ -3t  (ACn  + Cn  ) 
IF  q ' DPt  DINDt 


where  CDp  = 0.  01C335  is  the  overall  skin  friction  coefficient  (includes  wing, 
fuselage  and  tail). 

(viii)  Pitching  Moment  Coefficient 

The  pitching  moment  coefficient  about  the  wing  pivot  C ^ ^ (see  sketch  below) 
was  taken  as 


= 


cMp 


[ <ci/6'  > j - CL.  (bi/c)  sin^]  A Kt 


i=l 


+ CF  GLOVE  [xp/c-x  GLOVE/6] 

_ St  qt  fcr  cosQ;  +Cn  sino- 
S q c l ^ 1 


where 
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ACM6.  = -2(.03)(sine2-e2)sin5i2  - 2(.  O5)(sin0g  -93)sin6i3 

-2T1  (.17)(sin05  +05)sin6i5  - 2Tli6  (.  2O)(sin06  f 96>6iC 


cj/c  = 1.  0,  c2/c  = 0.  82,  bj/c  = 0.  525,  b2/c  - 1.  589, 
x p/c  = 0,  x GLOVE/c  = -0.  875,  h 0 = 2.  32,  St/S  = 0.  25,  qt/q  = 0.  9 


CDt  = 0.0055  + ACDpt  ' CDINDt 


(ix)  Buffet 


/pTP! 


= - (l-3ACn„  -3ACj 


TEi  10 


(x)  Boot  Bending  Moment 


'BB  At 


^ Z (bi/c)Ch.AKi 


Mbb  = 214,668  (MS)  CRp  q^  (inch  lbs) 

where  MS  is  the  model  scale,  (MS  = 1 for  the  full  scale  aircraft). 


5.2.3  Sample  Computations  with  PAM 

Sample  computations  were  made  with  PAM,  to  examine  the  realism  of  the  calculated 
non  linearities  and  of  the  individual  effects  of  the  various  actuators,  and  to  generate 
the  pitching  moment  tables  described  in  Section  5. 1.  3 for  making  simulation  runs.  The 
computations  included  the  parametric  variation  of  each  of  the  twelve  actuators  about 
the  nominal  value,  in  addition  to  tail  incidence  angle  and  angle  of  attack,  according 
to  the  following  table: 


NOMINAL  POSITION 


A11 

A12 

A13 

A14 

A15 

A16 

A21 

A22 

A23 

A24 

A25 

A26 

at 

500 

| 200 

200 

500 

500 

500 

333 

9 

PERTURBATIONS  FROM  NOMINAL 


AA-q  A Aj2  AA13  ^Ai4  AAjg  AAlf>  AA21  AA22  AA23  AA24  AA25  Aa26  Aat 


■ 

— 

-333  i 

-500 

-333 

| -333 

-3° 

B 

0 

i 

0 

1 

0 

. 

-1 

333 

333 

333 

0 

1 

( 667 

667  1 

400 

667 

667 

1 

j 

| 

i 

3 

Because  of  space  limitations,  only  a limited  number  of  parametric  variations 

listed  in  the  above  table  have  been  included  herein.  Representative  plots  showing  the 

effects  of  the  A,  . A,  A,  and  A,,  actuators  on  the  variations  of  C with  a , 

13  14  15  21  L 

Cq  with  Cl  . and  Cj^with  Cl  are  shown  in  Figures  5.3(a)  - 5.6(c).  Similar  parametric 

plots  for  the  A actuator  and  for  at  were  given  previously  in  Reference  21.  The 

shapes  of  the  nonlinearities  in  the  stall  region  shown  in  the  figures  appear  sufficiently 

similar  to  test  data  to  make  simulations  canned  out  with  PAM  a meaningful  test  of 

effects  of  nonlinearities  on  convergence. 

The  pitching  moment  tables  were  generated  from  the  same  nominal  actuator 

positions,  except  A and  A^  = 200  instead  of  5 00  counts  and  av  = 0°  instead  of  5°. 

14  24  1 

The  calculations  were  carried  out  with  tail  incidence  angles  Qj.  between  -24°  and  12°, 
and  for  angle  of  attack  a between  0°  and  16°.  The  I’esults  may  be  plotted  as  shown 
in  Figure  5.7,  from  w'hich  tabulated  pitching  moment  data  given  in  Table  5.  2 were 
prepared.  These  data  are  to  be  used  in  simulations  involving  inequality  constraints 
on  tail  trim  control,  as  discussed  in  Section  5. 1.3. 


58 


5.3 


COMPUTER  SIMULATIONS 


An  extensive  set  of  simulations  involving  the  optimization  computer  program  in 
conjunction  with  PAM  were  run  on  the  DEC  10  computer.  The  majority  of  the  runs, 
which  were  reported  originally  in  Reference  22  and  are  summarized  in  Tables  5.  3 and 
5.4,  dealt  with  effects  of  cycling,  testing  errors,  the  size  of  the  perturbations  in  the 
incremental  mode  (DA's),  etc. , on  convergence  for  problems  dealing  with  either 
minimizing  CD  or  maximizing  C^  A summary  of  the  "standard"  values  used  for  the 
components  of  DA,  DEV  and  TOL,  and  a list  of  the  "standard"  testing  errors  is  found 
in  Table  5.5.  The  inequality  constraints  were  on  trailing  edge  pressure  (buffet 
constraint)  and  on  differential  actuator  deflection  (twist  constraint)  as  follows: 

cPTEis  0,  Cp^^O 

| A 1 j “ A2 j | 5 200,  j =2,  3,  4,  5 

Equality  constraints  were  enforced  on  pitching  moment  (trim)  and  on  lift  coefficient, 
as  shown  in  Tables  5.3  and  5.4.  The  tail  incidence  angle  was  taken  as  a free 
independent  variable  to  be  optimized,  hence  the  constraint  on  C^  could  be  satisfied 
without  recourse  to  the  approximate  procedure  given  in  Section  5.  1.  3. 

Representative  results  from  these  runs  are  presented  in  Figures  5.  8 to  5. 12. 
Thus,  Figures  5.  8-5. 10  illustrate  the  effect  of  cycling  (averaging)  and  the  choice  of 
the  random  number  (seed)  for  generating  testing  errors  on  convergence  for  the  minimum 
Cq  problem.  Figure  5.8  depicts  the  minimum  value  of  CD  obtained  for  sets  of  runs 
with  1,  2 and  3 cycles.  Six  cases  are  shown  for  each  set,  corresponding  to  different 
seed  values  in  starting  the  random  number  generator.  Although  the  runs  with  standard 
testing  errors  consistently  gave  higher  values  of  C^MIN  than  those  without  testing 
errors,  it  was  somewhat  surprising  that  the  lowest  values  of  Cg>  we rc  found  for 
the  two  cycle  case,  rather  than  for  three  cycles,  as  was  expected.  This  result  is 
illustrated  in  Figure  5.9  which  shows  average  values  of  Cd  for  1,  2,  and  3 cycles 
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with  standard  testing  errors  and  without  testing  errors.  The  variations  of  Cl,  Cd  and 
CM  during  typical  runs  from  the  series  are  shown  in  Figure  5. 10. 

Figures  5. 11  and  5. 12  deal  with  maximizing  Cl-  All  constraints  are  the  same  as 
shown  in  Figure  5.8  , except  that  the  Cl  constraint  has  been  replaced  by  a constraint 
on  Co.  viz. , C[)  = 0. 11,  corresponding  to  a sustained  maneuver  condition  for  the 
F-lll  TACT  aircraft  at  M = 0. 90.  Figures. 11  shows  the  effect  of  cycling  on  convergence. 
The  maximum  CL  with  testing  errors  is  again  obtained  from  the  two-cycle  case,  but  is 
somewhat  below  that  found  without  testing  errors. 

The  study  shown  in  Figure  5.12  was  carried  out  to  investigate  effects  of  the 
size  of  testing  errors  and  of  the  size  of  the  actuator  perturbations  during  the  incremental 
mode  on  convergence.  Thus,  all  runs  were  made  with  twice  the  standard  testing  errors 
given  in  Table  5.5  For  the  1 cycle  case,  convergence  was  poor,  and  was  helped  some- 
what by  increasing  the  size  of  the  actuator  perturbations.  Increasing  the  number  of 
cycles  to  2,  3,  and  4 gave  similar  results  to  that  observed  previously. 

However,  the  uppermost  curve  in  Figure  5. 12  shows  that  good  convergence 
was  eventually  obtained  by  changing  the  number  of  cycles  and  actuator  perturbation 
sizes  in  combination.  This  was  accomplished  by  controlling  the  DEC-10  "interactively" 
from  a video  terminal  and  monitoring  the  resultant  convergence  properties.  This 
"interactive"  capability  has  been  added  to  the  optimization  program  and  will  also  be 
available  during  the  wind  tunnel  test,  because  of  the  similarity  between  the  hardware 
at  AEDC  and  at  GD /Electronics  Division  at  which  the  DEC-10  is  located.  A detailed 
discussion  of  the  "interactive"  operation  of  the  optimization  program  has  been  included 
in  Reference  4. 

Prior  to  the  transfer  of  the  computer  program  to  AEDC/ARO  (on  Dec.  1G,  1976), 
several  additional  simulations  were  performed  with  held  fixed  (frozen),  and  also 
with  other  actuators  frozen,  to  simulate  more  closely  the  testing  at  AEDC.  Run  215 
in  Figure  5.13  deals  with  minimizing  Cn-  The  equality  and  inequality  constraints, 
tolerances,  etc. , are  identical  to  those  for  the  previous  runs  with  at  variable, 
except  that  the  equality  constraint  CM  - 0 has  been  dropped  and  replaced  by  the 
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inequality  constraints  stated  by  Equations  5.22  and  5.23.  The  required  tabulated 
values  were  taken  from  Table  5.2. 

As  seen  from  Figure  5.13,  reasonable  convergence  was  demonstrated  with 
zero  testing  errors,  with  the  solution  giving  a considerable  lower  value  of 

0.  043)  than  had  been  obtained  for  the  trimmed  zero  testing  error  case  (CdmIN 
~ 0.049  from  Figures  5.8  and  5 9).  As  seen  also  in  Figure  5. 13,  the  pitching 
moment  coefficient  drifts  to  the  positive  maximum  limit  after  the  fifth  iteration,  and 
remains  at  this  bound  for  the  duration  of  the  run.  The  trace  of  Cm  during  the  run 
in  relation  to  the  bounding  values  from  the  tables  may  also  be  seen  from  Figure  5. 14. 

It  is  surprising  that  Cm  becomes  bounded  positively  rather  than  negatively,  since 
a positive  Cm  may  be  thought  of  as  adding  additional  drag.  This  anomaly  may  be 
due  to  the  position  of  the  C.G.  (assumed  at  the  pivot),  and  further  simulations  should 
be  carried  out  to  verify  the  cause  for  the  positive  Cm- 

Similar  optimization  runs  with  standard  testing  errors  have  prove  unsuccessful 
to  date.  The  runs  have  indicated  an  inability  to  satisfy  constraints.  A possible  cause 
of  the  difficulty  is  the  value  of  TOL  = . 003  used  in  satisfying  the  pitching  moment 
tables,  which  may  be  too  small  in  the  presence  of  testing  errors.  Additional  simu- 
lations and  fine  tuning  of  the  program  parameters  are  required  for  this  condition. 

Runs  involving  maximizing  Cl  with  c<t  frozen  and  with  pitching  moment 
tables  introduced  to  limit  the  maximum  and  minimum  Cm  values  were  also  attempted, 
but  were  unsuccessful.  These  runs  must  also  be  simulated  successfully  prior  to 
attempting  them  in  the  wind  tunnel. 

The  optimization  computer  program  was  modified  prior  to  shipment 
to  AEDC/ARO  for  up  to  a maximum  of  twenty  dependent  functions  and  constraints. 

It  was  also  demonstrated  in  a series  of  interactive  runs  that  various  actuators,  which 
might  not  be  present  in  the  full  scale  AFTI-111,  could  be  readily  locked  out.  These 
modifications  included  freezing  the  upper  surface  humping  and/or  leading  edge 
sharpness  systems,  not  deflecting  about  the  25%  chord  position,  and  testing  without 
differential  twist  along  the  leading  edge.  Although  the  optimization  program  may  be 
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used  with  C ^ (corrected  for  tail  trim  lift  according  to  Equations  5.25),  demonstrations 

of  the  equivalence  between  optimizing  on  CJ  (corrected)  with  the  tail  held  fixed,  and  on 

with  a self  trimming  tail,  were  not  obtained.  It  is  recommended  that  this 

equivalence  be  established  by  making  additional  simulations  prior  to  testing  at  AEDC 

» 

with  the  tabulated  pitching  moment  constraints  and  with  Cl  replaced  by  CL. 
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TABLE  5. 1 DICTIONARY  OF  COMMON  VARIABLES  IN  THE 
OPTIMIZATION  COMPUTER  PROGRAM 
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TABLE  5.  2 PITCHING  MOMENT  TABLES 

OBTAINED  wrm  PSEUDO-AERODYNAMIC 
MODULE  (PAM) 


ALP(I) 

CM  MIN  (I) 

CMNOM(I) 

CM  MAX  (I) 

_2 

-0. 176 

0.020 

0.310 

0 

-0. 140 

0.044 

0.336 

2 

-0. 120 

0.068 

0.364 

4 

-0.090 

0.092 

0.390 

6 

-0.060 

0.115 

0.418 

8 

-0.030 

0. 137 

0.450 

10 

-0.005 

0. 157 

0.480 

12 

+0. 025 

0.175 

0.510 

14 

0.070 

0. 188 

0.550 

16 

0. 190 

0.  244 

0.650 

LTAB  = 10 

a*  =0 
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200, 

200,  500,  333,  333, 

500,  200,  200, 

ij  = 11, 

12, 

13,  14,  15,  16, 

21,  22,  23, 

TABLE  5.3  RUN  LOG  - DEC-10  OPTIMIZATION  RUNS 


i.mna.inAiinQ  no  atqpXo  >1VA  A A c II  r / r.ll-Hl 

^ rw|Ofj  mnpinni  |o  ||VA  A A 7.  A A /'/•/'  A Cl  I -Oil 


TABLE  5.4  DEC -10  OPTIMIZATION 
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Figure  5. 1 Flowchart  of  wind  tunnel  operation  as  controlled  by 
the  optimization  computer  program 
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Figure  5.3a  Effect  of  Inboard  Leading  Edge  Deflection  on  Lift  Curve, 
Pseudo  Aerodynamic  Module 
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Figure  5. 3c  Effect  of  Inboard  Leading  Edge  Deflection  on  Pitching  Moment, 
Paeudo  Aerodynamic  Module 
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Figure  5.  4b  Effect  of  Inboard  Upper  Surface  Humping  on  Drag  Polar, 
Pseudo  Aerodynamic  Module 
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Figure  5.4c  Effect  of  Inboard  Upper  Surface  Humping  on  Pitching  Moment, 
Pseudo  Aerodynamic  Module 
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Figure  5.5b  Effect  of  Inboard  Trailing  Edge  Deflection  on  Drag  Polar, 
Pseudo  Aerodynamic  Module 
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Figure  5.  6a  Effect  of  Outboard  Leading  Edge  Radius  on  Lift, 
Pseudo  Aerodynamic  Module 
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Figure  5.  6c  Effect  of  Outboard  Leading  Edge  Radius  on  Pitching  Moment, 
Pseudo  Aerodynamic  Module 
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Figure  5.  9 Effect  of  Number  of  Cycles  on  Average  Value  of 
Minimum  CD  with  Standard  Testing  Errors 
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Figure  5. 13  Convergence  of  Drag  Minimization  with  Tabular  Pitching  Moment 
Constraints  During  Run  215 
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Figure  5. 14  Variation  of  Pitching  Moment  with  Relation  to 
Tabular  Rounds  During  Run  215 


R:n\  2_:  IT 


FR0L.E.Ki  ST  1E.KQ 
Cs-tfiM  TPCtE-S  FSO'.M  TAGLE_  -a 


l-in  IT  F-F.3/-I 


S-RuATlOU  J.  Z3 


C.M  VAX  I AT /OKI 

during  id  'll  S' 


(osgs) 


I O.lo 


6.0  AERODYNAMIC  ANALYSES 


The  purpose  of  the  aerodynamic  analyses  has  been  to  evaluate  the  relative 
correspondence  of  transonic  aerodynamic  characteristics  between  the  target 
airfoil  shapes,  to  which  the  SOFT  Wing  has  been  designed,  and  the  closest 
approximating  (actual)  shapes  achieved  with  the  model  wing.  The  target  shapes 
have  been  defined  at  the  0. 348  b/2  and  0.913  b/2  semispan  stations  and  in  the 
streamwise  direction  with  the  wing  swept  at  the  -A=  16°  position. 

The  list  of  target  and  actual  airfoils  from  among  which  the  comparisons 
will  be  made  have  been  listed  in  Table  6. 1.  The  target  airfoils  have  been  designated 
by  the  prefix  TACT  or  CSC,  whereas  the  actual  airfoils  have  been  designated  with 
the  prefix  SOFT.  The  ordinates  for  the  various  target  ai roils  in  the  local  chord 
plane  (untwisted)  are  listed  in  Tables  6.  2 and  6.  3 for  the  inboard  and  outboaixl  span 
stations,  respectively.  Section  profiles  have  been  plotted  in  Reference  4.  The 
corresponding  actual  airfoil  ordinates  must  awmit  the  completion  of  skin-on  calibra- 
tion of  the  wind  tunnel  model,  as  discussed  in  detail  in  Reference  4. 

The  transonic  airfoil  computations  have  been  carried  out  for  the  target 
airfoil  shapes  in  Tabic  6. 1 for  the  condition  M = 0.  80  and  Of  = 0°  (corresponding 
to  the  component  of  the  free  stream  normal  to  the  leading  edge  at  26°  and  M = 0.  89). 
An  available  F.  Bauer  version  of  the  Jameson  nonlinear  transonic  potential  theory 
program  was  used  (based  on  Reference  23  and  24).  The  program  was  supplied  to 
GDC  by  Mr.  Ray  Hicks  of  NASA  (ARL),  and  includes  boundary  layer  displacement 
effects  in  combination  with  relaxation  methods  to  solve  the  finite  difference  equations. 

To  save  on  computer  time,  and  because  only  a relative  comparison  of 
aerodynamic  characteristics  is  to  be  made,  only  the  first  iteration  in  the  relaxation 
procedure  was  taken  and  no  boundary  layer  displacement  effects  were  included. 

Figure  6. 1,  based  on  previous  inviscid  computations  which  were  made  with  the 
same  computer  program  for  the  airfoils  presented  in  Reference  10,  shows  that 
the  initial  iteration  approximates  lift  coefficient  to  within  6 C ^ ~ 0. 05,  at  least 
for  the  three  airfoils  presented. 
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The  single-iteration  inviscid  Bauer  program  computations  of  Cj  » C^  and 
for  the  target  airfoil  shapes  are  listed  in  Table  G.  4.  The  correspondence 
between  the  calculated  chordwise  pressure  distributions  for  the  TACT  W52  and 
creepless  supercritical  CSC  airfoil  sections  (airfoils  1 and  3 inboard,  and  2 and 
4 outboard)  has  been  indicated  in  Figures  6.  2 and  6.  3. 

Following  the  determination  of  the  actual  airfoil  shapes  it  is  planned  to 
perform  the  comparitive  computations  of  Cl*  Cd  , CM  and  Cp  distribution  versus 
chord.  A newer  faster  version  of  the  Bauer  program,  recently  received  from  NASA 
(ARL),  may  be  used  for  the  actual  airfoil  computations  (and  possibly  to  further 
evaluate  the  accuracy  of  the  single -iteration  inviscid  approximation).  However, 
the  "fast"  Bauer  program  has  as  yet  not  been  converted  to  run  at  GDC  on  the  CDC 
CYBER  70  (or  at  GDE  on  the  DEC  10). 


TABLE  6. 1 LIST  OF  AIRFOILS  FROM  AMONG  WHICH  COMPUTATIONS 
AND  COMPARISONS  WILL  BE  MADE  WITH  ONE-ITERATION 
IN  VIC  ID  BAUER  PROGRAM 


Status 

Case 

Mach 

Angle  of 

Inboard  Airfoil 

Outboard  Airfoil 

0- not  run 

1- run  completed 

No. 

No. 

Attack 

0.  348  b/2 

0. 913  b/2 

TACT,  W52 

CSC-  (t/c=0.  0958) 

SOFT- 1A  (Matches  #1) 

SO  FT- IB  (Matches  #3) 
TACT,  W52(TE  #5) 

TACT,  W52(TE  #6) 
TACT,  W52(6LE=5°) 
TACT,  W52(6lE=10°  ) 
TACT,  W52(5LE=15°  ) 


TACT,  W52 

CSC-(t/c=0.  06122) 

SOFT- 2A  (Matches  #2) 

SO  FT-  2 B (Matches  #4) 
TACT,  W52(TE  #5) 

TACT,  W52(TE  f! 6) 
TACT,  W52(6LE=5‘) 
TACT,  W52(6LE=10») 
TACT,  W52(  5 15°  ) 


SO  FT- 1C  (Matches  #9-1) 


SO FT-2C (Matches  #9-2) 


SOFT  lG(Matches  #9-9) 

| SOFT-2 G(Matches  4(9-10) 
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TABLE  6.4  SUMMARY  OF  TRANSONIC  FLOW  COMPUTATIONS 


One-Iteration,  Inviscid,  Bauer  Computation 
M=0.  80,  Qf=0* 


CL 

CD 

CM 

0. 9593 

0. 00534 

-0.2286 

0. 9994 

0. 00501 

-0.2385 

0. 9284 

0. 00421 

-0.2220 

0.  8478 

0. 00456 

-0. 1979 

1.4102 

0.01441 

-0.3243 

1.3852 

0. 01456 

-0.3231 

1.9690 

0. 03860 

-0.4438 

2. 0006 

0.04908 

-0.4636 

0. 9495 

0. 00281 

-0.2319 

0.9379 

0.00218 

-0. 2353 

0. 9275 

0.  00297 

-0.2384 

0.  9853 

0. 00043 

-0.2435 

0.  9673 

0. 00085 

-0.2490 

0.9498 

0. 00018 

-0.2527 

Airfoil 


TACT  (Inboard) 

TACT  putboard) 

CSC-  (Inboard) 

CSC-  (Outboard) 
TACT+TE  #5  (Inboard) 
TACT+TE#5  (Outboard) 
TACT+TE  #6  (Inboard) 
TACT+TE#6  (Outboard) 
TACT,  6lE=®°  (Inboard) 
TACT,  <5le  = 10°  (Inboard) 
TACT,3lE=15°  (Inboard) 

TACT,<5LE=5°  Putboard) 

TACT,(5le=10°  (Outboard) 
TACT, <5 LET  15°  (Outboard) 


Airfoils  from  Reference  10 
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Figure  6. 1 Variation  of  Lift  Coefficient  with  Iteration  Numhe 
for  Bauer  Computer  Program,  Inviscid  Analysis 


Symbol  Case  Airfoil  C^  Cj^  ^OL 

I TACT,  W52,  0. 348  b/2  0.9593  -0.2286  0.00534  -5.044 

II  CSC,  9.58%  0.9284  -0.222  0.00421  -4.878 


Figure  6. 2 Comparison  of  TACT  and  CSC  Airfoil  Characteristics  from 
Bauer  Version  of  Jameson  Program,  Inboard  Station 
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Figure  6.  3 Comparison  of  TACT  and  CSC  Airfoil  Characteristics  from 
Bauer  Version  of  Jameson  Program,  Outboard  Station 


7.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  present  interim  technical  report  has  dealt  with  progress  made  under  Phase  1 
of  the  Self  Optimizing  Flexible  Technology  Wing  wind  tunnel  program.  The  Phase  1 
investigation,  which  included  preliminary  and  detailed  design  of  a 1/G'th  scale  semi- 
span  model;  loads,  stress  and  flutter  analyses;  fabrication  of  a kinematic  model, 
the  computer  programming  of  an  optimization  program;  simulation  studies ; and 
various  aerodynamic  analyses  has  been  completed.  The  Phase  II  program,  which 
deals  with  model  fabrication  and  testing  at  AEDC,  is  now'  W'ell  underway,  anil  plans 
call  for  testing  in  mid  summer  of  1977.  Specific  conclusions  and  recommendations 
with  regard  to  the  Phase  I investigation  are  listed  below. 

(i)  Detailed  design  of  the  wind  tunnel  model  has  been  completed  and  fabrication 
is  currently  underway.  The  model  was  designed  with  the  planform,  thickness 
distribution  and  twist  distribution  of  the  TACT  I wing  and  includes  12  independent 
hydraulic  actuation  systems,  three  flexible  skin  panels  with  sliding  joints  and  a 
split  spar.  The  articulation  is  sufficiently  general  to  encompass  proposed  AFTI-111 
wing  contours,  as  well  as  transonic  wings  with  more  general  camber  and  thickness 
variations,  and  provides  for  conformal  spanwise  as  well  as  chordwise  changes 

in  wing  shape.  The  design  permits  increasing  the  number  of  independent 
actuators  to  20  and  inserting  small  chord  extensions  between  the  spars,  with  only 
minor  rework,  to  even  further  increase  the  generality  of  the  model. 

(ii)  All  design  objectives  appear  to  have  been  met,  except  that  problems  are 
anticipated  in  achieving  differential  twist  variations  along  the  span  (primarily 
along  the  leading  edge,  but  also  possibly  along  the  trailing  edge),  and  the 
articulation  may  be  severely  limited  for  the  inboard  trailing  actuators  when  the 
wing  is  swept  at  58°.  In  addition,  the  trailing  edge  thickness  was  compromised 
at  1%  chord  (instead  of  the  0.  75%  chord  desired),  and  fasteners  were  allowed  to 
protrude  above  the  skin,  because  of  problems  in  designing  the  sliding  trailing 
edge  joints 
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(iii)  The  stress  analysis  of  the  split  spar  design  showed  an  adequate  minimum 
safety  factor  SK-  5.  7,  based  on  ultimate,  for  the  primary  wing  structure  at  the 
nominal  test  dynamic  pressure  q-500  Ibs/ft^  Safety  factors  have  also  been 
shown  to  be  marginally  adequate  for  the  various  lever  arms  and  brackets  com- 
prising the  actuators  (lowest  SF-3.  9,  based  on  ultimate).  ()n  the  other  hand, 

it  has  proven  difficult  to  estimate  the  maximum  stress  levels  in  the  flexible 
skin  panels.  It  is  therefore  recommended  that  up  to  50  strain  gages  be 
implanted  on  the  skins  during  the  pretest  calibrations  and  the  maximum  skin 
deflection  limits  be  established. 

(iv)  The  wing  flutter  analysis  has  predicted  positive  aerodynamic  damping  at 

/ 2 

transonic  tost  conditions  for  q ^1250  lbs/ft  . Thus,  no  flutter  problems 
are  anticipated  for  the  actual  dynamic  pressure  levels(q~  500  lbs/ft“) 

However,  because  of  the  unusual  wing  design,  and  because  of  the  many 
simplifying  approximations  in  the  flutter  analysis,  it  is  recommended  that 
static  deflection  checks  be  made  during  the  pretest  calibrations  (to  verify 
the  stiffness  assumptions),  and  that  aerodynamic  damping  runs  at  gradually 
increasing  q levels  be  made  at  the  beginning  the  the  test. 

(v)  Preliminary  estimates  of  panel  flutter  have  shown  the  upper  rear  skin 
panel  to  be  somewhat  marginal.  It  is  therefore  recommended  that  two  strain 
gages  be  premanently  implanted  on  this  panel,  and  that  the  corresponding 
RMS  strain  reading  be  monitored  during  the  test  for  any  indication  of  panel 
flutter. 

(vi)  A wind  tunnel  scale  kinematic  model  with  10  actuators  has  been  fabricated 
and  was  used  to  verify  the  function  of  the  leading  and  trailing  edge  actuators  and 
the  form  of  the  wing  shapes  produced.  Conclusions  reached  through  operation 
of  the  kinematic  model  were  that  reasonably  good  transonic  (and  possible  low- 
speed)  airfoil  sections  could  be  produced  with  the  present  design,  and  that  the 
anticipated  level  of  articulation  should  be  obtainable,  except  that  wrinkling  of 
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the  leading  edge  skin  was  observed  when  differential  deflection  of  the  leading 
edge  actuators  was  attempted.  In  addition,  the  push  rod  (for  changing  leading 
edge  radius)  was  found  to  be  too  flexible,  and  was  also  observed  to  unport  from 
the  skin  when  the  leading  edge  was  deflected  downward.  The  design  of  the  push 
rod  has  been  modified  and  the  calibration  procedure  has  been  altered  to  accommo- 
date the  above  difficulties.  However,  problems  are  still  anticipated  for  leading 
edge  differential  twist. 

The  kinematic  model  may  also  be  used  for  estimating  actuator  force 
(and  hydraulic  pressure)  requirements  and  for  predicting  areas  of  skin  stress 
concentration  and  skin  stress  levels,  both  unloaded  and  loaded.  However, 
funding  and  scheduling  considerations  prevented  these  applications  from  being 
carried  out.  It  is  recommended  that  future  SOFT  Wing  programs  take  fuller 
advantage  of  the  analytical  tools  which  the  kinematic  model  offers. 

(vii)  The  2D  optimization  program  has  been  completely  rewritten  for  the  SOFT 
Wing  test.  The  new  optimization  program  allows  up  to  20  independent  vaiiables 
(actuators,  angle  of  attack,  tail  incidence)  to  be  involved  in  the  optimization, 
and  permits  the  merit  function  being  optimized  to  be  selected  from  among 
18  different  aerodynamic  and  geometric  functions  being  measured.  The 
remaining  17  functions  may  then  be  selected  as  equality  or  inequality  con- 
straints. The  program  features  a simple  means  of  freezing  (locking  out) 
various  independent  variables,  in  order  to  represent  wings  with  reduced 
articulation  and  includes  a provision  for  tabular  limits  on  pitching  moment 
when  optimizing  with  a fixed  tail.  The  program  may  be  controlled  inter- 
actively from  a remote  video  or  teletype  terminal . 

The  optimization  program  has  been  demonstrated  on  a DEC  10  computer, 
and  delivery  to  AEDC/ARO  has  been  carried  out. 


(viii)  A series  of  test  simulation  anti  check  out  runs  were  made  on  the  DEC  10 
by  linking  the  optimization  program  to  a Pseudo-Aerodynamic  Module  for 
computing  aerodynamic  coefficients.  The  majority  of  the  simulations  dealt  with 
maximizing  Cl  or  minimizing  Cj>  subject  to  equality  constraints  on  Cp  or  Cl> 
respectively,  and  on  Cm>  and  to  inequality  constraints  on  CpTE  ant' 

on  differential  actuator  deflection.  These  simulations,  which  were  made  with 
all  twelve  actuators  active  and  also  with  variable  tail  incidence  angle,  demon- 
strated reasonable  convergence  of  the  optimization  procedure,  even  with 
standard  testing  errors,  provided  that  the  various  control  parameters  could  be 
fine  timed  during  interactive  control. 

(ix)  Convergence  was  not  demonstrated  for  similar  optimization  problems 
with  testing  errors,  when  the  tail  was  held  fixed  and  the  equality  constraint 
on  Cm  was  replaced  by  tabular  inequality  constraints  to  limit  Cm  between 
values  for  which  adequate  trail  trim  control  was  available. 

(x)  It  is  recommended  that  the  simulations  be  completed  for  the  problems 
with  a fixed  tail  and  with  tabular  limits  on  Cm-  addition,  it  is  strongly 
recommended  that  no  optimization  problems  of  any  kind  be  attempted  in  the 
wind  tunnel,  until  adequate  simulation  and  fine  tuning  of  the  optimization 
parameters  has  been  achieved  including  effects  of  testing  errors  This  is 
required  to  insure  that  valuable  wind  tunnel  occupancy  time  will  not  be 
squandered  on  procedures  which  could  be  carried  out  at  very  low  cost  on 
the  DEC  10. 

(xi)  The  "fast''  Bauer  computer  program  should  be  converted  for  running 
on  the  DEC  10  and  used  to  evaluate  the  accuracy  of  the  one-iteration  inviscid 
approximation  made  in  the  aerodynamic  analyses. 
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